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Theory of Oscillation Type Viscometets: 
The Oscillating Cup’) 


ParT I 


By JosrerH Kestin and GorpDOoN FRANK NEWELL, Providence, R. I.?) 


Abstract 


This is the first of a series of papers dealing with the theory of oscillation 
type viscometers with particular emphasis on those in which a finite right 
circular cylinder (cup or disk) oscillates in contact with a fluid. The purpose is 
to obtain formulae which relate the density and viscosity of the fluid to the 
frequency and logarithmic decrement of the oscillation and the various physical 
properties of the suspension system. The study is to include an accurate analy- 
sis of the fluid motion in the vicinity of the edges since this is the main 
source of error in most existing theories. The first paper gives a general formu- 
lation of the problem and an exact solution (including transients) for the 
oscillating cup viscometer with the fluid inside the cup. 


List of Symbols 


The following is a partial list of symbols including only those that are used 
repeatedly: 


A surface of contact between fluid and cylinder; 
D(s) defined in equation (16) ; 
f(t) transient, equation (42) ; 

h half height of cylinder or height of liquid in an open cup; 
I moment of inertia of suspension system; 

I’ moment of inertia of fluid inside cup; 

I, modified Bessel function of order k; 

Jy Bessel function of order k; 

M(t) frictional moment exerted by the fluid; 

n normal to A outward from the fluid; 
A radial coordinate in the fluid; 


1) This work was supported by the United States Air Force through the Air Force Office of 
Scientific Research and Development Command under contract No. AF 18(600)1548. Part of the 
work was done while one of the authors (G. F. N.) was being supported by a grant from the Alfred 

- P, Sloan Foundation. ; 

-2) Brown University. 
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R radius of cylinder; 

Sm see equation (34) ; 

see equation (25); 

see equation (38); 

S, roots of equation (39) ; 

t time; 

T) natural period of oscillation ; 

w dimensionless form of 2, see equation (12) ; 
2 vertical coordinate in fluid; 

a(t) angular displacement of body; 

a initial angular displacement of body; 
6 boundary layer thickness (v/a )"?; 

A logarithmic decrement ; 

A, logarithmic decrement in vacuum; 

7 dimensionless vertical coordinate 2/6; 
N) dimensionless height h/d; 

fs zeros of Jy, Jy(t4;) = 9; 

» kinematic viscosity of fluid; 

Eg dimensionless radial coordinate 7/6; 

&, dimensionless radius of cylinder R/0d; 
ig density of fluid; 

03 dimensionless time WT; 

o circular frequency of oscillation ; 

@, circular frequency in vacuum; 

92 angular velocity of the fluid; 

* complex conjugate; 

bar over symbol denotes the Laplace transform of the function. 


1. Introduction 


The use of oscillating systems for the measurement of the viscosity of fluids 
leads to great simplicity of design and a high accuracy of measurement. The 
only measurements required are those of length, mass, and time all of which 
can be made very accurately, the errors being reducible to 1 part per 1000. : 

The types of systems in most common use consist of the following arrange- 
ment: A cylindrically symmetric body is suspended by a torsion wire along 
its axis so as to have a certain observable natural period and logarithmic decre- 
ment in a vacuum. If the body is then placed in contact with a fluid, the 
oscillation of the body will also impart motion to the fluid with two principal 
consequences. (1) The fluid is viscous and some of the energy is dissipated in 
heat. This gives rise to an increased value of the logarithmic decrement of 
oscillation for the combined system. (2) The fluid has mass and an effective 
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moment of inertia which causes the system as.a whole to have a longer period 
of oscillation. 

The changes in frequency and decrement caused by the fluid depend only 
upon the density and viscosity of the fluid and measurable properties of the 
suspension system and, therefore, in principle, one could infer the viscosity and 
density of the fluid from measurements of the frequency and decrement. In 
most cases the density can be measured directly with panthe accuracy so that 
only the viscosity is of interest. 

The main difficulty that has prevented the use of this type of system for 
precision viscosity measurements is the absence of sufficiently accurate theore- 
tical relations between the observed frequency and/or decrement of oscillation 
_and the unknown viscosity. Typical applications using existing approximate 
formulae have resulted in errors ranging from a few percent to as much as 
100 % in the calculated viscosity. 

The following is the first in a series of papers in which the theory of some 
of the more common types of oscillating systems will be extended so as to give, 
wherever possible, formulae that are as accurate as is necessary in order to 
make full use of the accuracy of existing data. 

The specific types of oscillating systems that have been proposed or are in — 
use are classified as follows: 

(1) The solid body is a disk or circular cylinder of finite height suspended in 
a virtually infinite fluid [1-3]8). 

(2) An arrangement as in (1), except that the solid body is a sphere [1, 4, 5]. 

(3) The solid is a disk again but the disk oscillates in a plane parallel with 
two fixed plates placed above and below the disk so as to retard the motion of 
the fluid situated in the region between the plates [3]. 

(4) The body is a hollow cylinder or sphere and the fluid is contained inside 
the body. The region outside the body may be a vacuum or it may be a fluid. 
In the latter case the fluid motion outside can be treated separately and is 
equivalent to one of the above cases [5—7]. 

The literature dealing with either experimental or theoretical aspects of 
these problems is quite lengthy. The references noted above are some of the 
more recent publications most of which contain more extensive bibliographies 
of earlier work. 

Of the cases listed above, only the spherical bodies lead to relatively 
simple exact formulae. These have already been treated elsewhere [5] and will 
not be considered here. Unfortunately there are some experimental difficulties 
involved in the design and manufacture of this type of viscometer and it has 
not been used as much as the others. 

The remaining cases all have one principal difficulty i in common. None of 
_ the existing theories except [6] accurately account for the peculiar behavior 


3) Numbers in brackets refer to References, page 465. 
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of the fluid in the vicinity of the edges where the vertical and horizontal sur- 
faces of the cylinder meet. Most of the following work will be directed toward 
the analysis of these edge effects which are the principal source of error in 
present theories. 

The first case to be treated in this and in the second paper is the hollow 
cylinder with the fluid inside and a vacuum outside: the oscillating cup visco- 
meter. This will be studied in considerable detail mainly because this is the 
only one of the cylindrical shaped bodies that has been treated exactly. The 
analysis of this serves a dual purpose of being an end in itself and also as a 
testing ground for various approximation methods. An exact solution for the 
cup has actually been obtained and used by SHVIDKOVSKII ef al. [6| but they 
give only one form of the solution which in some situations is very difficult to use. — 

The following sections also give a general formulation of the mathematical 
problems common to all the other viscometers to be studied later. 


2. Formulation of the Problem 


The formulation and notation used in this section follows closely that 
described in more detail in [8, 1, 3, 4]. 

The motion of the body is described by the equation 

2 

Tog [HO 4 2A, “+ (1 + AB) afe)] = MO), (1) 
in which «(z) is the angular displacement of the body from equilibrium, 
Wy = 2n/Ty is the natural angular frequency of oscillation in a vacuum, T, is 
the corresponding period of oscillation, t = gt is a dimensionless unit of time ?, 
I is the moment of inertia of the solid body, J w@ is therefore the spring 
constant, M(z) is the frictional moment exerted by the fluid, and A, denotes 
the logarithmic decrement in a vacuum [for M(t) = 0]. 

The analysis of equation (1) is facilitated by the use of the Laplace trans- 
form. We assume that initially the system is displaced from its equilibrium and 
then released so that ; 
da(t) 


a(t)'= a, and 7 


=i) aiore ye == (0) (2) 


If we denote the Laplace transform of any function f(t) by 


co 


f(s) = [eo He) ae, (3) 


0 


then «(s) satisfies the equation 


[(s “te dae oF 1] «(s) = 
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In order to determine M(r), or its transform M(s), we must investigate the 
motion of the fluid, and for this we make the usual assumption that no secon- 
dary motion is developed, i. e., we omit the non-linear terms in the Navier- 
Stokes equations. The present analysis is, strictly speaking, applicable to zero 
amplitude oscillations or to data extrapolated to zero amplitude. In practice, 
however, the amplitude is usually small enough so that no extrapolation is 
necessary. 


Figure 1 


System of coordinates. 


The fluid motion can be described in terms of its angular velocity Q which 
in a polar-cylindrical system of coordinates 7, », z (Figure 1), satisfies the dif- 
ferential equation [5], 


(5) 


where v is the kinematic viscosity of the fluid. The boundary conditions to be 
satisfied by 2 are 


BO hg {07 302 020 
Ot ~ ? Or? a in Oe le 


Qy, z,t) = a at the boundary, r= R, z=0 and z= 2h, | 
(6) 

OGs, Oy =O. 

The first condition implies that there is no slip at the boundary between the 

fluid and the oscillating body. The second condition implies that the fluid is 


initially at rest. 
The torque M(r) is related to the fluid motion through the equation 


M(z) = —0 Bhd Selo (7) 
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in which @ is the density of the fluid, A denotes the surface of contact between 
the fluid and the body, 7 is the normal to A from the fluid into the body and 


dois an element of surface area. 
It is convenient to introduce dimensionless space coordinates defined by 


R h ie 


along with the dimensionless time t = @ ¢. The length 6 is, in some sense, an 
average boundary layer thickness. 
Equation (5) now becomes 


I0% 2020 tes due 00 lieeED 
dee OEY? FOE oOneae (9) 


If there is no fluid motion initially, the Laplace transform of equation (9) gives 


— 1020 3- 02 022 


$0 = of EOE On? (10) 
and equation (6) gives the boundary condition 
O = Wy (S&H —%) onA. (11) 
We can also make Q non-dimensional by choosing | 
v= ff ; (12) 


Wy (S % — H&) 


The fluid motion is then described by the completely dimensionless equation 


0?w 3 0? 
Ssw= ja tert One? w= w(E, 7, s) (13) 
with the boundary condition 
We. 4; Ss) = ona se (14) | 


ae : 
M (s) can be expressed in terms of w and the result substituted in equation . 
(4) to give 
a(s)) ot ieee 


% Ss  sf(s +4)? +14 D(s))’ (15) 


Ds) = +2 Teale ga eg (16) 


in which ” and do are also measured in units of é. 


with 
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By inversion of the Laplace transform, a(t) can be written as the integral 


— _%_ st {1 we Ai 
a(7) hes ¥: Ace Ames TO oF 


along any vertical contour C in the right-hand half of the complex plane. 

The mathematical problem is to solve equations (13) and (14) for w, or, 
strictly speaking, for just dw/dn at the surface only, to determine D(s) from 
equation (16) and then to substitute in equation (17) to find a(t). 

The goal is primarily that of finding accurate formulae for «(t) which are 
not so cumbersome as to be useless. It should be kept in mind that the ulti- 
mate purpose is to infer the value of viscosity from the observed behavior of 
a(t). Unless a(t) is relatively simple or at least approximately so, the problem 
of finding the viscosity by fitting theoretical curves for various values of » 
with experimentally observed curves is impractical. 

Present experimental techniques for using such an instrument are based 
upon the assumption that «(z) will be very nearly a damped oscillation of the 
type 

a(t) ~ e—4°* cos(wt + y) (18) 
and the values of viscosity are to be inferred from measurements of w or A 
or both. 


3. Exact Solution for the Cup 


In treating the oscillating cup, there are two slightly different physical 
situations to be considered. In the first scheme the cup is closed at the top 
with the aid of a lid which remains in contact with the fluid (e. g., a gas). In 
the second, the fluid (e. g., a liquid) possesses a free surface. In the first case 
the boundary conditions are as given in equation (14), whereas in the second 
case we have 


= 0 de er cen iia. (19) 
where 7, denotes the dimensionless height of the column of liquid, on condition 
that surface tension is neglected. ’ 

These two cases are mathematically identical, because in the first case (cup 
of height 2 7,), the plane at 7 = % is a plane of symmetry along which condi- 
tion (19) is satisfied. It is, therefore, sufficient to solve the problem for the 
- lower half of the cup only; the drag on the total surface will be just double that 
on the bottom half. ) 

The simplest procedure for solving equations like equation (13) are to use 
separation of variables or to expand in some suitable complete set of orthogonal 
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functions. Actually these two basic procedures are quite similar because the 
most suitable set of functions in which to make expansions are those which 
naturally arise from the separation of variables. There are, however, several 
modifications of these basic schemes which lead to several different forms of 
the solution. Since the methods in question are fairly well known, we will 
simply outline the principal steps. Five different forms of the solution each of 
which will later be particularly useful in some special situations were obtained 
by the five methods described below. 
(i) We notice that the function 


w'(é, ; s) _ w(é, UB s) pd 
satisfies the inhomogeneous differential equation 


0?w’ Ow 07w’ Ds 
TL re Meee a. On? Sw =1§ (20) 


with homogeneous boundary conditions 
w= (0 on.A - (21) 


The set of functions €-! J,(u,; €/&), with J, denoting the Bessel function of 
first order, satisfy the differential equation 


ha Wee AUR ACESS 
age te Hillier age ae (22) 


which arises from equation (20) with s = 0 by separation of variables. If we 
define the parameters yu; so that these functions will vanish at = &,, then 


Silus) = 9, 


i.e., the y;’s are the zeros of J;. These functions form a complete orthogonal 
set and any function which is continuous in the interval 0 < &é < €, and van- 
ishes at € = & can be expanded in a convergent series of such functions. In 
particular we can write 


(ee) 


w'(E, ns) = 372, s) EB (EE). 3) 


j=l 


If we substitute this into equation (20) and make a similar expansion of the 
right-hand side of equation (20), we obtain an equation which can be satisfied 
only if the 0; are solutions of the ordinary differential equations 


hat Za Be abe SL ’ 
a | bi( a) ia Mh; A) (24) 
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> 
ue ey tS. 


(25) 
Since w’ vanishes on the boundary we must require that 


b(0, s) = b;(2 yo, s) =0. (26) 
The solution of equations (24) and (26) is 


My St Jo(us) 


' —2&s 
6; (n, s) = .; 
and so from equation (23) we obtain 


cosh sy (49 — ae 1] 
cosh sy, 19 


ha 1 Jf coshs, (7) — 7) 
w é; ? S$} = 1 lp 2 Ss a 5 | f a ae , 
(5, , s) a ica 


eNO cS) Sis E/E) 
cosh sy %o 1] Faia lie (27) 


If this is substituted into equations (16) and integrated term by term, one 
obtains the final result, 


{, _ Ss [, — tanhsy mo ]\ 
arateaeal Ip 


7 Su No (28) 
I’ represents the moment of inertia of the fluid inside the cup; 
V' = 5 m0 m & 55 = — nok Rt (29) 
for a fluid with a free surface at 7 = 7) and 
IPE=TO Ro, 0 =—woh _ (29b) 
for a full cup with a lid at 7 = 2 %p. 


It is also possible to obtain an alternative form of equation (29) by using 
the identity 8 D’ w;* = 1. This gives 


nb ie th tanhs, 7 
ero 3 poesia AMAIA 28 
vis § Z aa 3 = bi se i Su No | ee 


(ii) If we had a very long cylinder, we know that, for most values of 7, 
w(é, 7, s) is nearly independent of 7. We would then have 


0?w 3 Ow : 
. —_ ~ ie 1 
OF se ee | a al 0 with w(&>, 4, s) 
The solution of this equation is 


wee celal 


ET, (Vs &) ’ 
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in which J, is the modified Bessel function of order 1. This is, of course, the 
exact solution for a cylinder of infinite height. 
These facts suggest that we try to find the difference function 
M é, I, (Vs €) 


yg! = — — 


EI, (/s &o) 


which satisfies the same differential equation as w but the boundary conditions 


Uf = 1 = N a Eo Ty (V/s E) 
w" (Eo, 9, s)=90, w"(&, 0,s) = w"(E, 2%, s) =1 EL, (Vs &) (30) 


As in (i), we again seek to find an expansion of w” of the type 
= HA My &g 1 
w"(&, 7, 8) Tei ies ahintee” &o Ns (31) 


Substitution of this into the differential equation gives a set of homogeneous 
equations for the 6”, namely 


a? 2) ys 
gat ~ Sep se D. 
We require that 0;'(0, s) = 0;(2 mo, s) so that 6/(y, s) must be of the form 


" ey cosh Su (n + No) 
b; (yn, s) = 8; (0, s) cosh a aa | (32) 
The values of the coefficients 6/'(0, s) are finally determined by substituting 
equation (32) into equation (31) and forcing w” to satisfy the conditions given 
in equation (30). Integration of the solution gives the result 


- 4 I, (V/s &>) 8s tanh,(su 7) 
DI S\ce=S= 0 (Su No 
y Vs oh, (Vs) 0 2 Hy Sit |. (33) 


This expression is exactly equivalent to equation (28) and could have been 
derived from equation (28) if one knew the appropriate expansion identities 
to transform from one to the other. The method and the result given by 
equation (33) are essentially the same as given in [6]. 

(ii) We can apply a procedure which differs from that given in (i) only in 
that the role of the € and 7 coordinates are interchanged. The set of functions 
sin[(2 m+ 1) 7/29] bear the same relation to equation (20) as did the 
functions €-1 J(u; &/&). They are solutions of the equation 


(a ee ce 
dn? Bair ae 7m ie 


oe 
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which arises from equation (20) by separation of variables in the same way as 
did equation (22). 

These functions vanish for 7 = 0 and have a vanishing derivative at nH = No, 
provided m is an integer. They form a complete set and any continuous function 
satisfying the same boundary conditions can be expanded in a convergent 
series (Fourier series) of these functions. 

Corresponding to equation (23) we seek a solution in the form 


co 


w'(E,n, 8) = 3" ai,(é, s) sin CMF 27 


m=0 2 No 


In analogy with equations (24), (25) we obtain 


ee es 4s 
Vat Fae 5 ASE Salen car esr 
with 
Cire. ¥ (2m+1)2 2 
f= s+ ( se ) (34) 


and the a;, must satisfy the boundary conditions 
(E>. 8)= 0, -@,(0, s) finite. 


The solutions of these equations finally lead to the result analogous to equation 


(28), ee 
ecxien ee: 1 im pins | 
55 ae wit 2, (2m +1)? : Sm  Ly(Sm &o) J’ >) 
and a summation of the leading terms gives the alternative form 
= I’ 3 8 r + T3(Sng 7) 
Dis) = S457 8° a (mei 8, Lise (35a) 


(iv) The following method bears the same relation to (ii) as (iii) did to (i). 
If we have a very squat cylinder, the € dependenice of w is only slight for most 
values of . We expect that 


gr 8 ~0 with w(&,0) = w(f,2) =1, 


wo ~ £OSBVS (to = 7) 
cosh//s 19 


the solution for two infinite plates. 
This suggests that we try to evaluate 
Mw — cosh)'s (10 =) ® 


he cosh//s No 
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We expand this in a Fourier series and determine the coefficients in a manner 
analogous to that of (ii). The only difference is that the 7 and & coordinates 
exchange roles in the same manner as in (i) and (iii). The procedure is straight- 
forward and gives the result, 


I’ | tanh (//s x 32s . 1 Llemaeolal 
D(s) = 3? = t ne us idee aE, ¢ ya eee vj: (36) 
(v) The final method is, in essence, a combination of (i) and (ii). The set 
of product functions & 1 J,(u; &/&) sin[(2 m+ 1) 2/2] form a complete 
orthogonal set for expansions of any function of the two variables € and 7 that 
vanishes on the boundary of the cup. The product functions are, of course, the 
‘elementary’ solutions of the differential equation obtained by separation of 
variables. 
If we make a double series expansion of w”’ 


w'(é, n, S mea o> Gin é-1 mA ae ) ain (Z ear! 1) : 


substitute this into equation (20) and make a similar expansion of the right- 
hand side of equation (20), we obtain immediately the coefficients 4;,, which 
define the solution. Actually one can proceed in a number of other ways to the 
same result, namely 


BIE wee i 
OD ee rae Tn? ras pe (2 m 4)? (Ss — S5mq) oo 
with 
eke oflem ly 2 
Sim = ( 2 | Bre ee 


We now have at our disposal a sizeable collection of exact forms for D(s), 
namely equations (28), (28a), (33), (35), (35a), (36), and (37), however, any 
comparison of usefulness of the various forms for D(s) or any approximations 
must be considered in conjunction with the evaluation of the integral in equa- 
tion (17) for the actual motion of the cup, «(t). 

Formally, at least, «(t) can be evaluated by residue theory since the only 
singularities of the integrand are poles. This gives 


ee (1 + Aj) exp(S, 7) 
at) mo 2 ae oy 


in which S; are the roots of the equation 
(Sz + Ap)? + 1+ D(S;) = 0, (40) 


i, e., the poles of the integrand. D’(S,) denotes the derivative of D(s) at's= oe 
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and the k-th term of the series represents the residue of the integrand at the 
point s\=.S,. 

If I’/I = 0, then D(s) = 0 and equation (40) reduces to a simple quadratic 
with only two solutions; 


(Sy + AP 1 =O; “S,= += 4,, (41a) 
and 
a(t) = e~4°* &, (cost + Ay sint) . (41b) 
This, of course, represents the motion of the cup in the absence of any fluid. 
If we now consider what happens for 0 < I'/I <1, we see that regardless 
of how small /’/I may be, D(s) cannot be small for all values of s because it 
has infinitely many poles at the points s;,, [equations (37), (38)]. Near each 
of these poles which lie on the negative real axis there is one point where D(s) 
will assume any given value of order 1 or larger and in particular there is one 
point S; near each s;,, where D(S;,) has a value that satisfies equation (40). 
Since all quantities in equation (40) are real for real values of S;, it is easy to 
prove that each of the new roots S, is also a negative real number. 
Besides the infinitely many negative roots which lie near the s;,, we still 
have two roots which for J’/J<1 lie near the points + 7 — A). Since the 
complex conjugate (denoted by *) of equation (40) gives 


(S® 4Ag)? 1 DiS?) = 05 


it follows that if S;, is a root of equation (40), so also is S¥. Any root of equation 
(40) must, therefore, either be real or one of a complex conjugate pair and the 
two roots near + 7 — A, must still form a complex conjugate pair even for 
Pj + 0. 

Despite the fact that J’ >0 causes a singular behavior in the roots of 
equation (40), «(z) approaches equation (41b) continuously as I’ > 0. It is 
easy to see that the terms of equation (39) associated with the two complex 
roots vary continuously into those for J’ = 0 and give equation (41b), so it is 
only necessary to see that as J’ > 0, the contribution from the infinitely many 
negative roots contribute nothing to «(t). This is true because at each negative 
root S,, D(S;,) is at least of order 1 and since S;, is very near a pole of D(s), 
D'(S;,) is very large, of order I/I’>1 to be exact. D’(S;,) dominates the denomi- 
nator of equation (39) and causes each of the new terms to be of order J Ait 
For I’/I <1, «(t), therefore, is of the form 


a(t) = e7~4°* a (A coswt + Bsinw t) + f(t) (42) 


with /(t) of order I'/J and the various other parameters differing from those of 
“equation (41b) by quantities of order I’/J. 
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The A and w in equation (42) are obtained from the two complex roots S;, 
and S* by writing them as 


S,=o@(-4 +2). (43) 


Since equation (40) is complex, one may choose to consider the two real 
equations obtained by taking the real and imaginary parts of equation (40), 


namely 
Re D [w (-A + i)] = —1+ w? — (Aw — A)? (44a) 


Im D [ow (—A + 2)] = £2 w (Aw — Ap) . (44-b) 


Many of the results obtained above for I’/I <1 are modified only slightly 
for larger values of J’/I. The transcendental equation for the S;, 1s analytic 
in I’/I and the values of S, will therefore also be analytic functions of, Fa, 
The roots must still be either real or form complex conjugate pairs and those 
which were real for I’/I <1 will be real for all values of I’/I. 

To investigate the real roots further, we consider the behavior of the ‘unc- 
tion D(s)/s? along the real axis. We observe from equations (37) and (38) that: 
(a) D(s)/s? has poles at sj, and the s;,, are independent of I'/I, (b) the residues 
at these poles are all positive, (c) the derivative of D(s)/s? is everywhere nega- 
tive, and (d) D(s)/s2->0 for s+ co. Condition (b) implies that D(s)/s* 
becomes positively infinite as one approaches any pole from the right and 
negatively infinite from the left. This along with (c) implies that between every 
two consecutive poles on the real line, D(s)/s? is monotone decreasing from 
-+oo to —oo and (d) guarantees that D(s)/s? > 0 to the right of the largest s,,,. 

The S, are solutions’ of the equation ; 


D(Sy) = oft Ay 2 
Sta, eae (ee 


both sides of which are sketched in Figure 2 for real S,. Figure 2 is not drawn 
to scale but serves only to illustrate the shape of the curves as described above. 
The desired values of S;, are at the intersections of the two curves. 

It is apparent from the graph or the above description of its properties, that - 
there is a value of S; between every pair of consecutive poles of D(s). As I’/I 
increases, the values of S, move continuously to the left but always remain 
trapped between the values of s;,,. 

The two roots which were at ++ « — Ay for I'/I = 0 also move continuously 
as I’/I increases and it is clear that «(t) can always be formally separated into 
two parts as in equation (42). The part f(t) which we call the transient part 
consists of an infinite series of terms each of which has a time dependence 
given by the factor exp(S;,t) with S, < 0. The magnitude of the largest S;, 
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Figure 2 


Graphical scheme for locating the roots S,. 


(smallest | S,|), S,, is a measure of the rate at which the slowest transient term 
decays. Figure 2 shows that S, is less than the largest s;,, S91, but larger than 
the second largest s;,,, which may be either spp or s,,, thus, 


a Ee | - ce = So2 ' 3 2 
ise ye inne ers Moanin Ieee (45) 
12 thy ae = 


4. Survey of Practical Applications 


The motion of the cup is described by equation (42) and to guarantee the 
possibility of accurate experimental measurements of w and A and accurate | 
calculation of y and o from them, we must insist that the following conditions 
be satisfied. (a) The transient /(t) in equation (42) is either uniformly small 

compared with %», decays rapidly as compared with the decay of the main 
oscillation or varies so slowly compared with the main oscillation that it is 
essentially a constant. In other words, we want f(t) to be such that its presence 
does not obstruct the observation of the main oscillation. (b) A must be small 


448 JosepuH Kerstin and GORDON FRANK NEWELL ZAMP 


enough so that one can observe at least a few oscillations before the amplitude 
becomes too small to measure. As a practical matter this means A < 1/10 
approximately. (c) We must have usable formulae from which » and @ can be 
calculated from A and w to an accuracy consistent with the accuracy of the 
experimental data. (d) A and w must be sufficiently sensitive to the values of 
y and 9 so that small errors in the observation will not translate into large 
errors in the calculated values of v. In addition to these, there are other 
practical limitations imposed by the availability of proper materials, instruments 
etc., but these problems will not be considered here. 

All the above conditions relate in one way or another to the properties of 
the function D(s) which is also a function of the parameters I’/I, ny and 5. I’/J, 
equation (29), depends upon the density of the gas and on the suspension 
system but not on v; 7%) and &), equation (8), depend upon » and the geometry 
but not on 9; whereas 7)/£) = h/R is independent of the fluid. A determination 
of I’/I and 7» or &, from observed values of A and using equation (44) is thus 
equivalent to a determination of o and », respectively. To test whether or not 
in any particular situation the above conditions are satisfied, one must have 
at least some preliminary estimate of 9 and v and thus also I’/I, &) and %p. 

D(s) s-2 I/I’ is a function only of s & and s 7 so that in essence I’/I scales 
the magnitude of D(s), & or 7% scale its argument whereas 79/& = h/R deter- 
mine the functional form. All properties of the cup are therefore naturally 
classified according to the magnitude of I’/J, & or m7) and h/R, thus according 
to the inertia of the fluid, the size of the cup and its shape. Actually a more 
suitable measure of ‘size’ is the quantity | so,| = [2/2 7]? + [u,/&>]? and we 
shall designate a cup as ‘large’ or ‘small’ as this quantity is small or large 
compared with 1. A small cup is one for which the boundary layer thickness 
is large compared with either R or / and it gives rise to the physical situation 
in which the fluid moves with the cup almost as a rigid body. A large cup is 
one for which the boundary layer thickness is small compared with both R and 
hand gives rise to a situation in which only a small layer of fluid with a distance 
of order 6 from the walls of the cup is set in motion. 

The various consequences of having large, small or intermediate values of 
I’/I, h/R and ‘size’ will be considered individually in Part II where suitable 


approximate forms for D(s) will be obtained for each case along with precise . 


conditions under which the conditions above listed are satisfied, particularly 
(c) and (d). We can, however, make certain preliminary estimates here based 
upon the qualitative properties of D(s). 

If I’/I <1, as would probably be true if the fluid were a gas, we see from 


the discussion of equation (42) that we have the generally favorable conditions — 


that f(t) = O(e% I'/I) <a, A — Ay = O(I'/I) , and w = 1+ O(I"/I). This, for 


sufficiently small J’/I, automatically guarantees that conditions (a) and (b) are | 


satisfied regardless of the shape or size of the cup. 


EE ef 
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We see from equation (45) that the slowest transient is proportional to 
exp(Sp t) < exp (so, t) and decreases by a factor of order exp(—2 | Sox|) in 
a time corresponding to one natural period of oscillation. Since we describe 
a cup as small if |so,|>1, it follows that for the small cup the transient 
decays very rapidly, in a time much less than one natural period and this will 
guarantee that condition (a) is satisfied regardless of the size of J’/I or of 
h/R. Actually the condition |s9,| S>1 is much stronger than necessary because 
if condition (b) is satisfied, it would follow that the decay time of the main 
oscillation is longer than the actual period of oscillation which in turn is 
always larger than the natural period. For the transient to decay rapidly 
compared with the main oscillation, it would, as a practical matter, suffice if 
the transient decayed to 1/e of the value in one natural period or less. Thus we 
would consider condition (a) to be satisfied if |s9,| > 1/22. 

By arguments similar to the above, we would conclude that for a large cup, 
the slowest transient decays very slowly, but this, in itself, is of little importance. 
One must consider all the transient terms together particularly those with a 
shorter decay time. The analysis of the large cup requires special treatment in 
almost all aspects and will be discussed in Part II. We will find, however, that 
condition (a) usually is satisfied for the large cup by virtue of having both a 
slowly decaying and small amplitude transient. The one situation in which 
condition (a) may give trouble is that in which I’/J is of order one or larger 
and the dimensions of the cup are such that the decay time of the transient is 
comparable with the natural period (|s9,| equal to about 1/10). Actually in 
this situation where condition (a) is likely to be violated, we will find that 
most of the other conditions give trouble also. 


Zusammenfassung 


Diese Abhandlung ist der erste Bericht aus einer Veréffentlichungsreihe, die 
sich mit der Theorie schwingender Viskosimeter befasst. Dabei sind besonders 
die behandelt, bei denen ein endlicher Kreiszylinder (Becher oder Scheibe) von 
einem fliissigen Medium gefiillt oder umgeben schwingt. Die Aufgabe der Unter- 
suchung ist es, die Abhangigkeit der Frequenz und des logarithmischen Dekre- 
mentes der Schwingung von der Dichte und der Viskositat der Fliissigkeit und 
von den verschiedenen physikalischen Gréssen des Aufhangungssystems zu 
finden. Weiterhin ist eine genaue Analyse der Fliissigkeitsbewegung in den Rand- 
zonen eingeschlossen, da diese die Hauptursache der Fehler in den meisten der vor- 
handenen Theorien war. Die erste Abhandlung gibt die grundlegende Formulie- 
rung des Problemes und eine exakte Lésung (einschliesslich der Einschwingzeit) 
fiir das Viskosimeter, bei dem ein mit der zu messenden Flissigkeit gefiillter 
schwingender Becher verwendet wird. 


(Received: June 17, 1957.) 
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Theory of Oscillation Type Viscometerts : 
The Oscillating Cup 
Part II 


By Davip Avan BeckwitH!) and Gorpon FRANK NEWELL’), Providence, R. I. 


Abstract 


The following is a continuation of Part I. The emphasis in this part is on 
inquiring how the frequency and decrement of the oscillation are related to the 
viscosity and density of the fluid for cups of various shapes and sizes. In par- 
ticular, estimates are made of the errors that would result from calculations 
of the viscosity and density from observed values of the frequency and. decre- 
ment and known experimental uncertainties of the latter. Methods are also 
described which enable one to perform these calculations explicitly. 


1. Introduction 


In Part I) an exact solution of the oscillating cup was obtained and some 
of the qualitative features of this solution, particularly the transient part, 
were discussed. The main objective now is to analyze in more detail the proper- 
ties of the main oscillation and to determine, among other things, if conditions 
(b), (c) and (d) of section I,4. can be satisfied. Although exact solutions are 
available, they are difficult to use because v and 0 are related to A and w 
through a pair of simultaneous transcendental equations, equation (1,44), 
involving D[w(— A +. 2)] which has been found only in the form of infinite 
series. The solution of these equations, therefore, requires special consideration. 

The properties of the oscillating cup depend upon its size, shape, and the’ 
moment of inertia of the fluid: Because various limiting cases show quite 
different properties and require different mathematical treatments, it is most 
convenient to consider these limiting cases separately. We consider first 
various shapes of small cups for which | s,| > 1, equation (I,45) ; then the large 

1) Brown University, now at Bell Telephone Laboratories, Murray Hill, N. J. 

2) Brown University. 


*) This is a continuation of Part I, preceding paper, and will employ the same notation. 
Equation numbers, figures, etc. referring to Part I will be designated by a prefix I. 
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cups, |Sox| <1; then the class of systems with I'/J <1; and finally any 
remaining cases. 

Much of the following analysis will center around the problem of finding 
accurate but manageable approximations for D(s) when s = w (+. i — A). In 
the absence of any fluid motion, @ = 1 and A = Ag, we anticipate that under 
most conditions in which we will be interested, we can assume at least that 
Ay <1, A <1, o is of order 1 and |s| = w|+i— A| is therefore also of 
order 1. D(s) s~* J/I’ is a function only of s & and s 72 and it will usually be 
implied without further comment that |s &| and |s | are of the same order 
of magnitude as & and 7% respectively for s = w (+ 7 — A). 


2. The Small Cup 


By ‘small’ cup, we mean 


rol= [a] (eT 0 


which implies that either 7 or &) or both are small compared with 1. We have 
already seen in section I,4. that this gives a rapidly decaying transient (even 
for |so,| 21/22) so that the transient needs no further comment and we 
turn our attention to finding convenient expressions for D(s) for |s| = O(1). 

We note that equation (1) implies also that |s;m| = |Sq,:|=> |s| for all 
j’s and m’s with |s| of order 1. An expansion of equation (1,37) in powers of s 
would actually be a series in s/s;,, for various values of 7 and m and should 
therefore converge rapidly. It is also a power series in y~? because 


vom SOREN Le (SMT LS @ 


is proportional to y for all values of 7 and m. This expansion of equation 
(1,37) gives ; 


with 


Regardless of the values of R and h, the first term (j= 1, m= 0) of equation 
(4) gives at least the correct order of magnitude of a, and therefore successive 
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terms in equation (3) differ approximately by a factor s/so, = O(1/s9). The 
fractional error in | D(s)| resulting from terminating the series, equation (3), 
after J terms is roughly |1/s,,|!+1. If we want this fractional error to be of 
order 10-? or less, for example, and are willing to compute no more than five 
terms in equation (3), then we must require that 


| Sox| Py 3 (5) 


which we accept as the practical limitation of usefulness of the small cup 
approximation replacing equation (1). 

We have obtained a rapidly convergent series in equation (3) at the expense 
of creating a doubly infinite series representation of the coefficients a, in 
equation (4). One should notice, however, that the a, depend only upon / and 
R, in fact a, can be written as h' or R! times a function of only h/R. This is of 
considerable practical importance because one would use the same cup to 
measure the viscosity of many fluids under various conditions of pressure, 
temperature, etc. The a, do not depend upon the fluid and need be evaluated 
only once for any given cup. 

If the given values of R and h are of comparable size, the series in equation 
(4) converges rapidly enough so that a numerical term by term evaluation 
presents no practical difficulty. Furthermore the rate of convergence increases 
rapidly with increasing /. If R/A is very large or very small, a term by term 
evaluation of equation (4) is rather tedious, however, and we shall seek better 
ways of evaluating the a, One could obtain suitable approximations directly 
from equation (4) but only by using considerable caution because one sees very 


soon that a convergent power series expansion of a, in powers of R/A or — 


h/R does not exist. A much better scheme of finding the a, for R/h <1 or 
k/h> Lis to use one of the other forms of D(s) derived in section I,3. 

We consider first the case R/h <1. One would expect that equation (1,28) 
or, preferably, equation (1,33) would be a good starting point as suggested 
partly by the manner in which these equations were derived. The expansion 
difficulties mentioned above can be attributed to the function tanhs,, 79 of 
eae (1,28) or (1,33). From the definition of S,, equation (1,25), we see 
that 


2 2 5 hX2 : 
Sig = (8&5 + 43) (Fe) S> 1 for all 7, (6) 
and therefore 
tanhs, 79 = 1—2exp(—2s, No) tree. (7) 


The term 2exp(—2s, 7) has approximately the order of magnitude 
exp(—,;h/R) S exp(—y,h/R) and clearly does not have a power series 
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expansion in R/h but it is asymptotically smaller than any finite power of 
k/h. By treating separately these exponentially small terms, one can expand 
the remaining part of D(s), equation (I,33) is powers of s and R/h thereby 
obtaining for a, the expressions 


m(s Sart) + (8) Sea gt 


i = +S [1 +0 (=)]. 


The quantities of the type 3’ u;” can be evaluated from tables of the j;. One 
finds that 


ee IO AR hes 101 4 10S 
(8a) 
ees eG) Suet 07 x10. 


The expressions for the a, in equation (8) are not exact but the exponentially 
small term in a, is an approximation that gives at least the order of magnitude 
of the most important error term. The fractional error in evaluating v will be 
essentially equal to the fractional error in a,. If we neglect this exponential 
term, the relative error in a, will be approximately (A/h) (0-3) exp (—3-8 h/R). 
This serves mainly to define the range of applicability of equation (8). Although 
we have assumed R/h <1, we see that even for R/h < 1/2, the relative error 
in a, is less than 10~4 and for R/A < 2/3, the relative error is less than 10-3 
which is about as accurate as presently needed for most measurements. 

The other a,;, / > 1 in equation (8) also contain exponentially small errors 
of similar type but, except for these, the expected series expaiuisions of the a, 
in powers of R/h conveniently degenerate into just two terms. 

To evaluate the a, when h/R <1, we use a scheme analogous to that describ- 
ed above for R/h <1. This time we start with equation (I, 35) or (I,36) and 
note that | s,, &)|S> 1 for all m, equation (I,34). We must, therefore, use an 
asymptotic expansion of the modified Bessel functions as the analogue of 
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equation (7). This leads to an asymptotic expansion of the a,, namely 


peo ee eye 


2h4 h 6 (16)? h4 [,. (z) 4 

17h® hh 30 (16)? A8 (F) oo | 
Katee TLte Bis) Paylin [és SE a | 

62 hs 1382 hi0 

A — ee = 
4 735 |: 8 550957 
in which 
oe . 
=v (2m4+1, (9a) 
m=0 


t,=1-:0045, ¢,=1-:0015, ¢,=1-0005, ¢,=1-0001. 


The accuracy of equation (9) is inherently limited by the asymptotic nature 
of the expansions of the a, but as a practical matter, it is more seriously limited 
by the use of only a few terms. If we want to compute a, to a relative accuracy 
of 10-3, then we must require /R <1/5, arather strong restriction on geometry. 

Figure 3 shows the approximate regions of the &, 7) plane in which the 
various formulae apply. The small cup region is where equation (3) applies and 
cases (a), (b) and (c) denote respectively the regions in which the a, are most 
easily evaluated using equation (8), (9) and (4), respectively. 

The physical significance of the first few terms of equation (3) is quite clear. 
For a small cup, the fluid moves with the cylinder almost as a rigid body and 
consequently the main effect of the fluid is to change the moment of inertia 
of the system from J to J + I’. If, indeed, we substitute only the first term 
of equation (3) into (1,40) or (1,44) and for simplicity neglect A? terms, we 
find, as expected, that 


i 1/2 
on + = (A ~ A, w) (10a) 
and so 
Ji 1 — w? 
Pile: Ot, et (10b) 


The viscosity does not enter the picture unless we take at least two terms of 
equation (3). To the next approximation we find 


Oy | 2 eT) (A Ay) 
v VET GD 


(11a) 
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Figure 3 


Separation of the 1), plane into regions in which various formulae of the text apply. Broken 
line, A may: is the approximate location of the maximum of A with respect to v. 


or ; 


A~ (4, +5 


Wo ie i, 1/2 
a eww lee ac ; te) 


The second term of equation (3) has relatively little effect on the value of w; 
the main effect of the viscosity is to produce an increase in the decrement due 
to the fact that the nearly periodic motion of the fluid is slightly out of phase 
with that of the cup. 

Additional terms in equation (3) can be attributed to the fact that the 
effective moment of inertia of the fluid is slightly less than J’ if there is a 
finite boundary layer thickness. 

Equations (10) and (11) permit us to test whether or not the small cup 
satisfies conditions (b), (c) and (d) listed at the beginning of section I,4. 
Condition (b), that A < 1/10, is easily satisfied because in equation (11b), 
a,(@,/v) is of order 1/| so, | which we already have assumed to be less than 
about 1/3, equation (5). The product J’ I”/?/(J + I’)? has an absolute maximum 
value of 2/3 /3 and J¥/?/(I + I’)? <1, therefore A < A, + 1/15 for all values 
of I’ and J. To guarantee that 4 < 1/10, it is only necessary that A) be suffi- 
ciently small. 
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Whether or not conditions (c) and (d) can be satisfied depends to some 
extent on how we choose to calculate » and @ from A and w. We anticipate 
the possibility that @, and therefore also J’, can be measured by other means 
and that y is actually the only unknown. Potentially a measurement of either w 
or A would suffice to determine » if is already known. A calculation of » from 
w and J’ would, however, necessitate measuring the small viscosity dependent 
corrections to equation (10) and consequently the inferred value of » would be 
very sensitive to the slightest errors in w and J’. If, on the other hand, we 
should calculate » from measured values of A and I’, we see from equations 
(10a) and (11a) that v1 is nearly proportional to A — A w and, if 4) @ <A, 
the fractional error in the computed value of v would be about equal to the 
fractional error in A (assuming that A is the least accurate measurement). If 
we should choose to determine both A and w from measurements of both A 
and w then we expect the errors in vy and @ to be comparable with the errors in 
A— A, and w?— 1 respectively. Since both decrements and changes in 
frequency can be measured quite accurately, the small cup is potentially a 
suitable instrument for determining both » and @. 

The final problem now is to devise ways of actually calculating y (and @) 
from A (and w) using equations (I,44a), (I1,44b) and (3). We will not give any 
more explicit formulas than these because such would only make the problem 
appear more difficult than it really is. One can, however, solve these equations 
explicitely for vy and J’ in terms of A and . IJ’ is eliminated by taking the ratio 
of equations (1,44a) and (I,44b). The resulting transcendental equation for y 
can then be solved in series. The fact that y-1is nearly proportional to (A — A)o) 
means that the series expansion for D(s) in powers of y~1 can be inverted to 
give y-! as a power series in (4 — w A,). Having determined y—1, one then 
substitutes this into equation (I,44a) and evaluates J’ as a similar series. 

In view of the fact that such formulae will be rather awkward, it may be 
simpler to solve the original equations numerically by a method of successive 
approximations starting with a trial solution determined by equations (10b) 
and (11a). 

If the value of I’ is known from other sources more accurately than it can 
be computed from w and JA, one should either eliminate w from equations 
(I,44a) and (I,44b) or solve them simultaneously for w and y instead of I’ | 
and y. In any case the amount of work required depends strongly on the rate 
of convergence of equation (3), the accuracy desired, etc. Space does not permit 
a detailed discussion of the numerical solution for all possible situations. 


3. The Large Cup 


The large cup is defined by |so,| <1 which obtains only if both Ny > 1 
and &) > 1. This suggests that one should try to obtain a power series expansion 
of D(s) in powers of yo * and & +. Since D(s) depends upon No and €, only in the 
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combinations s 75 = s h? w/v and s & = s R® we/v, such a power series would in 
effect be a power series in » or s~!. Unfortunately any attempt to make a 
Taylor series expansion of D(s) in integer powers of s~! fails almost immediately. 

We will sketch two different ways of finding an asymptotic form for D(s). 
The first method is based upon using the exact expression for D(s) given by 
equation (I,36); the second method is based upon an asymptotic solution of 
the original differential equation for w, equation (1,13). In both cases we will 
compute only enough terms to indicate the procedure and simply quote the result 
obtained by extension of the same methods. For a more detailed account see [9]. 

Starting from equation (I,36), we make a large argument expansion of 


tanh (Vs Mo) and J; (S &>). The lowest approximation gives 


D(s)~s? rig OES be aay 1 i (12) 


T \ Yin,” he oy Om +s 


From equation (1,34), we see that s>,~s at least for m< IVs nol. By substi- 
tuting, s,,° = s~*? + (s,,3 — s~3) and noting that D7’ (2 m +1)-? = 22/8, we 
obtain the equivalent expression 


Disjw 52 2 | af 


2 ca et 
| (12a) 


4 3 

a ent eet 

The infinite series represents a small contribution as compared with the first 
two terms but the evaluation of it and similar types of terms in higher order 
approximations require special care. For m< \Vs |, the numerators in the 
series are small but also nearly proportional to (2 m + 1)2/ng. The (2 m + 1)? 
cancels the one in denominator giving terms that are nearly independent of m. 
For all m, the terms of the series are slowly varying functions of m and to sum 
the series we must use the Euler-MacLaurin summation formula which yields 


is oe 1 ‘a AX US ee ay SORE 2 EL (13) 
a (2 m + 1)? sé 410 ; LN (see Be 2V5 no ; 
This is correct to any finite power of 79 1 but there is an error proportional to a 
negative exponential in 7. 
Substitution of equation (13) into (12a) gives three terms of the asymptotic 
expansion. A similar treatment of the higher order terms gives the result 


ibe + 6 3 3 
2 — : [= ate 
D(s) =S vi | sll? Ey s & ale 2 3/2 &3 = Dy 52 & i} 
i 16 ) 8 
Pains Pes 14 
T Sg wsiony | SB ng m S*E} No Oe 


= 2 rig? s exp (2.51 my) +--+}. 


458 Davin ALAN BeckwitH and GoRDON FRANK NEWELL ZAMP 


The first line is a series in powers of (s' &)-1 whereas the second line is 
(s'? m»)-1 times a similar series. There are no terms containing 79°, 7 = 2 but 
there are exponentially small terms the most important of which is contained 
in the third line. The ‘large cup’ region of Figure 3 shows the range of € and 
no where equation (14) gives D(s) to within about 0-1% for |s| of order 1. 

The second method of deriving equation (14) shows the physical significance 
of the various terms. It is clear from the physical meaning of a ‘boundary layer 
thickness’ that for £1 and 7,1, only a relatively small layer of fluid 
(thickness of order 1 <&p, 7 in units of 6) near the surfaces of the cup take 
part in the motion. Except within a distance of order 1 from the edge where 
the surfaces meet at right angles, the motion of the fluid near the verticle 
walls of the cup is only slightly influenced by the motion of the fluid near the 
bottom (or top) surfaces of the cup and vice versa. It is convenient, therefore, 
to consider the fluid motion as a linear superposition of three flows (a) the 
flow that would exist near the walls of a cylinder of infinite height (with no 
bottom), (6) the flow that would exist near a plane of infinite radius (with no 
cylinder walls), and (c) a flow localized near the intersection of the cylinder 
and the plane. Since D(s) is a linear function of w, D(s) is the sum of the 
contributions from each of these flows. 

The solutions for (a) and (6) are contained in section I,3. and these contri- 
bute to D(s) the two terms 


D(s) ~s? ge | 41, (V/s Eo) tanh (//s 7) (15) 


: Weep, (V/s &) Vs No 


An asymptotic expansion of the first term yields all the terms in the first line 
of equation (14) which we can, therefore, interpret as due to the drag of the 
fluid on the vertical wall. The first of these terms differs from the leading term 
of the small cup expansion, equation (3), by a factor of order 4 >! which one 
can explain by the fact that a moving layer of fluid of unit thickness near the 
wall has a moment of inertia equal to about 4 &1J’ as compared with the 
moment of inertia I’ of the entire fluid which is set in motion when we have 
a small cup. 

A large argument expansion of the second term of equation (15) gives the . 
first term of the second line of equation (14) which we interpret as due to the 
drag of the fluid on the flat surfaces and it gives also the exponentially small 
term of the third line of equation (14) which one may think of as due to an 
interference between the flows near the top and bottom surfaces. The factor 
no I’ that appears in these terms represents the moment of inertia of the 
fluid layer set in motion by the flat surfaces. 

The remaining terms in the second line of equation (14) must be due to the 


edge correction. The largest of these terms indeed contains a factor I’ ae ie 
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which has the same order of magnitude as the moment of inertia of the fluid 
within unit distance from the edge. These terms can be obtained by considering 
the fluid flow in the immediate vicinity of the edge. If one denotes by w, the 
difference (c) between the total flow w and the two flows (a) and (b) described 
above, then w satisfies the same differential equation as w but different bounda- 
ry conditions. An asymptotic solution for w, can be derived directly from the 
differential equation by making use of the fact that w, is small except near the 
intersection of a plane and a cyliader of relatively small curvature. This proce- 
dure is straight-forward but too long to describe in detail here [9]. 

We now consider the various requirements listed in section 1,4. for the 
cup to be suitable as a viscometer. We shall find it necessary, among other 


things, to postulate aoe ‘ 
ri eg)< v6 


So 


The analysis is greatly simplified by anticipating this result and making use of it. 

The arguments presented in section I, 4. regarding the suitability of a visco- 
meter when J’/J <1 apply also here if only equation (16) is true. This arises 
because D(s) is actually of order (I'/Z) (4 &1-+ 7’) for a large cup instead of 
order J’/I as assumed in section I,4. or in other terms because I’ (4 1 + 79 ') 
is the ‘effective’ moment of inertia of the fluid rather than I’. 

From this we conclude that the transient f(r), equation (1,42), is of order 
(I’/T) (4 & ++ 57). It is also possible to find simple approximations for /(r) 
using equation (16). One can do this by substituting the asymptotic form for 
D(s) directly into equation (I, 17), and deforming the contour of integration 
to the left of the two complex roots that give the main oscillation so as to lie 
along both sides of the negative real axis. This can be done provided t < & 
and 7% despite the fact that the asymptotic approximation for D(s) is incorrect 
on the negative real axis. A second method is to use equation (38) but to 
integrate over the densely distributed roots S; on the negative real axis. The 
result of this procedure is [9] 


[o,e) 


—% I’ (4 ty ge ae het ce ands ae ome Le 
f(t) ~ T [eet | = / (y2 +1)? or t< 5 and 19 ( ) 
and from this 
ut (A 1 . . ; 
Kz (= ie 7) 2a vl ns oe oe oe tele 


provided equation (16) is valid. 

This shows that the transient f(r) is always negative (relative to %») and has 
a monotone decreasing amplitude. For later times, t of order &2 or 7, T is of 
the same order of magnitude as the decay time of the slowest transient. The 
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transient is then dominated by the exponential decay of the slowest transient 
just as was the case with the small cup on a much shorter time scale. _ 

If equation (16) is not valid the transient will have a relatively large 
amplitude and such a complicated time dependence that it would be very 
difficult to observe the main oscillation. Equation (16) is therefore a necessary 
and sufficient condition for condition (a) of section I,4. to be satisfied. 

If we use equation (16) again and the consequence of it that s ~7 for the 
main oscillation, we find by substituting equation (14) into equation (I,44a 
and b) that to a lowest approximation 


il TRS iach 1 1 I’ (4 1 y \1/2 
l1-o~A eae nel I le ctere) 2y2 - (Gectesplilae) vl 

From this we see that equation (16) is also necessary and sufficient for 
condition (b) of section I,4. (namely 4 <1) to be satisfied provided that A, 
is also small. Equation (18) also indicates the accuracy with which various 
quantities can be calculated and the most important aspect to observe is that 
1 — w and A — A, are nearly equal. 

Unless one can measure both 1— w and A — A, so accurately that their 
difference can be inferred with significant accuracy (an unlikely possibility), 
an observation of both 1— qm and A — A, is redundant. An observation of 
either would to a first approximation be a measure only of the product J’ »¥?. 
If, however, o@ and thus J’ can be measured by other means with sufficient 
accuracy, then the fractional error in the computed value of » would be twice 
that of A — A, or of 1 — w, whichever quantity is used by itself as a basis 
for computing v. Potentially, both 1 — w and A — A, can be measured accu- 
rately, but in any particular situation one should use only that measurement 
which is the more accurate. ; 

Again we shall not try to give explicit formulas for y in terms of I’ and 1 — w 
or v in terms of I’ and A — A, to an accuracy consistent with the accuracy of 
equation (14). The procedure for doing this and the results are awkward but 
elementary. The series expansion for D(s), effectively in powers of v/?, can be 
inverted to give v? as a power series either in A — A, or in 1 — a. If one is to 
measure A and J’, one should eliminate w from equations (44a, b) whereas if 
one measures A and J’, one should eliminate A from equations (44a and b). 
Even though one may be able conveniently to measure both w — 1 and A — Ag 
if 1 — m can be measured more accurately than A — A), the theory implies 
that 4 — A, can be computed from 1 — w and J’ more accurately than it can 


be measured and conversely, if A — A, can be measured more accurately than 
1—o. 


4. Intermediate Size Cup 


In the previous two sections, we have considered all cases except those for 
which the smaller of €) and 7 is of order 1, the unshaded area of Figure 3. 
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The analysis of the remaining cases is complicated by the fact that theie is 
no suitable power series expansion of D(s) in powers of y~! or y"?, We shall find, 


however, that the behavior of the ‘intermediate size’ cup is very complicated 
unless 


7<! (19) 


and the analysis of this section will be based upon the consequences of this 
assumption. Some of the methods to be used here can also be used further to 
simplify the analysis of the small and large cups as well when equation (19) holds. 
We have already seen in section 1,4. that if equation (19) is valid then f(z), 
1 — w and A — Aj are all of order J’/J. The transient should give no serious 
trouble as a consequence of the fact that it has a small amplitude and will 
decay at a moderately fast rate for |s | 2 1/2. We will be more concerned 
with how 1 — w and A — A, depend upon the properties of the fluid. Because 
of the greater complexity of the present problem, we shall not try to describe 
it in quite as much detail as the previous cases but will simply describe the 
qualitative properties and indicate how one would proceed to an accurate 
analysis. 
. The starting point of this discussion is the pair of equations (I, 44a, b). For 
small J’/J, A and w can be expanded in a Taylor series in powers of J’/I. The 
leading term of this expansion gives 


@—1=4 ReD (41) + ol(7)'] (20a) 


A — Ay = +4 ImD (41) + ol(7)] (20b) 


The complete series is easily obtained and will involve various powers of 
D(- 1), which is of order J’/I, and various derivatives of D(s) at s = + 2. The 
simplifying feature of equation (20) is that it involves D(s) only at s = 2. 
D(-+ 1) is, however, still a function of & and my) and thus of the unknown vy and 
the known values of h and R. 
We obtain directly from equation (1,37), the expressions 


ReD (+1) = ——} +2 64? x wy? (2m +1)7 (1+ s;,)~" (21a) 

rag ee 7 64 a? a if @m+1) s+ Sie (21b) 

+ImD (-47) = Soa hg” (Ao YE Se) (Eb Soltis eters) 
q,m 


in which s;,, is proportional to », equation (2). 
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Each of the individual terms of equation (21b) is a monotone decreasing 
function of » and so, therefore, is the sum. Thus to order (J’/I)* at least, 
1 — w is a monotone decreasing function of » with 1 — w > 0 for y > 0 (limit 
of large cup) and 1 — w > (I'/I)/2 for vy > 0 (limit of small cup). For any given 
values of hand R, 1 — wis simply I'/J times a function of y only. There is nothing 
very unusual about this function and if the given values of K and h are of 
comparable size, it could be tabulated as a function of v quite easily from a 
term by term evaluation of the series. 

Each of the terms of equation (21c) vanishes for » > 0 (s;,, > 0) and for 
y > 00 (S;m —> oo). The (7, m) term has a single maximum with respect to v at 
= —1 thus at the point 


pose (Ee a oat (22 


Sjm 


The larger 7 and m, the smaller the value of y at which the 7, m term has its — 
maximum. Thus to order (I’/I)?, at least, A — Ay vanishes for y > 0 and for 
y oo as we already observed in the previous sections. It is also clear, however, 
that A — A, has a single maximum with respect to v and it occurs at a slightly 
lower value of y than that of equation (22) with 7 = 1, m = 0. If we denote by 
Yman the value of »y where A — A, is a maximum with respect to y for fixed R, 
hand I’, then 


2 ay eo | 
Pnas S09 \ eee + Ey (23) 


This is not only an upper bound on the value of y,,q,; it is also a fairly close © 
estimate. For R/h > 0, equation (23) implies that R(@9/?maz)!? > fy ~ 3-8 as 
compared with a correct value of 4-3 computed by Priss [6]. The broken line — 
of Figure 3 shows, approximately, the curve in the &, 7) plane where A has a 
maximum with respect to v for fixed R, # and I’. 

We can also obtain from equation (21c) simple upper and lower bounds on ~ 
the maximum value of A, Aimax. Amaz is larger than the maximum value of — 
the first term alone but smaller than the sum of the terms evaluated at ee | 
respective maximum values. This gives 


a a Oa ee 
qieh the IE Nip * 


eS (24) 
An exact calculation of Ajnq, for h/R >0 gives Amar — A, ~ I'/5:81 [6], 
which is about midway between the limits of equation (24). Fouagen (24) holds 
for arbitrary values of h and R, however: 
This is all the information we need to test the suitability of a cup for 
measurements of viscosity. Equation (24) shows that we can guarantee 


, 
tain icine a tl ct ee 
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A S1/10 [condition (b) of section I,4.] for arbitrary geometry if only 
1 
ae (5) 


and Ay is sufficiently small. This defines more precisely the range implied by 
equation (19). 

In regard to sensitivity of J’ and y to measurements of w and A, the existence 
of a maximum of J as a function of » poses new problems. If we should try to 
calculate » from measured values of J’ and A (but not w), we would find that » 
was not only a double valued function of A but also rather insensitive to the 
value of A if y should lie near »,,,,. One must, therefore, be careful to avoid 
trying to compute » from I’ and J if vis near y,,4,. Except for this region of y, the 
fractional errors in y will be at least comparable with (but probably somewhat 

aed than) the fractional errors in A — A, and I’. 

A calculation of y from measured values of J’ and 1 — w presents no special 
problems for the intermediate size cup, in fact, v has its greatest sensitivity to 
@ in this range. One should note, however, that fractional errors in y are 
comparable with fractional errors in 1 — w, not o itself. 

A calculation of both y and I’ from A and @ presents no problems, in prin- 
ciple. In effect, the value of J’ is inferred mainly from the value of A which is 
nearly proportional to Z’/J and v from the value of (1 — w)/(4 — Ay) which is 
nearly independent of J’/J. In this scheme, the fractional error in I’ is compa- 
rable with that of A — A, whereas the fractional error in » is comparable with 
that of A — A, or 1 — w whichever is the larger. 

The final problem of actually calculating y explicitly from any two of the 
quantities J’, A and w is somewhat more difficult than in the cases of the small 

_ or large cup but it is still not beyond reasonable limits of complexity. 

If we were to use a method of successive approximations (equivalent to a 
series of corrections in successive powers of J'/J) the first approximation would 
be the actual solution of equations (20a, b). If we measure w and J’ (but not 4), 
we eliminate A by considering only equation (20a) and solving it for v; if we 
measure A and J’ (but not w) we eliminate w by using only equation (20b) ; and 
if we measure w and A (but not J’) we eliminate J’, to order (I’/I)?, by dividing 
equations (20a, b). D(-+ 7) I/I’ depends only upon k, A and » and for any 
given cup (fixed R and h) it, therefore, depends only upon ». 

From a graph or tables of Re D (+ 1) I/I’, Im D(+ 4) I/I' and their ratio 
as a function of v, or better yet the inverse of these functions, the value of » can 
be found immediately to order I'/I. 

The tabulation of D(-- 2) can be done using equations (21a, b, c) if R and 

h are comparable size. These series will converge, in this case, fast enough to 
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be practical particularly in view of the fact that the various terms of each 
serics al! have the same functional form and involve only very simple operations. 

If R/hS> 1, it would be easier to compute the real and imaginary parts of 
D(-+- 7) from equation (1,33). The real and imaginary parts of the Bessel 
functions with arguments proportional to Vei can be expressed in terms of 
the tabulated functions ber and bei whereas the function tanh(s, 7%) in 
equation (1,33) can be represented by its large argument expansion since 
Su No ~ fy h/R> 1. If on the other hand h/R <1, it would be easier to find 
the real and imaginary parts of D(-+ 7) using equation (I, 36) and the asymptotic 

expansion of the Bessel functions for s,, By ~(2m+1)a R/2h>1. 

Successive corrections of y to higher orders in I’/J can be obtained from a 
series expansion in powers of I'/I, 1 — w or 1 — A depending upon which of 
these are ‘knowns’ and which are eliminated from the equations. Additional 
tables for derivatives of D(s) at s = +7 might be helpful but not necessary 
for these corrections. 

A detailed treatment of the various mathematical steps particularly for 
R/h <1 is described by Priss [6] although it is possible to organize the calcu- 
lations somewhat more efficiently than described there. 


5. Conclusions 


We have described here all the possible situations in which a cup filled with 
a fluid might be used successfully as a viscometer. Based upon the criteria 
listed in section I,4., we have found virtually no restrictions on the use of a_ 
‘small cup’, equation (1). For the large cup we found that equation (16) was 
a necessary condition for usefulness but that one could not easily determine 
both vy and @ from A and w. For the intermediate size cup we found equation — 
(25) was sufficient to satisfy most requirements although one cannot always 
compute » accurately from observations of only A and I’. Equation (25) is 
also a necessary condition for the intermediate size cup to be useful because if 
this condition is violated the transient. will be of relatively large amplitude, 
the decrement will be too large and the theoretical calculations exceptionally 
tedious, thus it violates virtually all of the conditions for suitability. 

Because of the many ways in which one might actually perform the final 
calculations of » (and perhaps @) from observations of A and/or w, we have 
avoided giving details of this. We have, however, in each case indicated some 
methods which would not involve an unreasonable amount of work. 

In view of the great flexibility of design and the potentially high accuracy 
that can be achieved, it appears that an oscillating cup should compete favorably 
with other types of viscometers particularly for measurements of the viscosity 
of liquids. Although some use has already been made of such an instrument, 
its potentialities far exceed its present applications. 


a re art 
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Zusammenfassung 


Diese Notiz ist eine Fortsetzung des ersten Teiles. Hierin wird insbesondere 
die Frage behandelt, wie die Frequenz und das Dekrement der Schwingung von 
der Viskositat und Dichte der Fliissigkeit, bei Verwendung von Bechern ver- 
schiedener Gestalt und Grésse, abhangt. Vor allem wurden die Fehler abge- 
schatzt, die bei der Berechnung der Viskositaét und der Dichte aus den beobach- 
_ teten Frequenz- und Dampfungswerten auftreten kénnen und wenn letztere mit 
bekannter experimenteller Ungenauigkeit bestimmt wurde. Weiterhin sind 
Methoden angegeben, wie diese Berechnung leicht erfiillt werden kann. 
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Stresses in a Plate Under Tension and Containing an Infinite 
Row of Semi-Circular Notches 


By Axrra Atsumt, Sendai, Japan’) 


Introduction 


The purpose of this paper is to determine the stresses in a semi-infinite 
plate with a row of an infinite number of semi-circular notches when it receives 
a uniform tension JT. The problem of determining the stresses in the neigh- 
bourhood of a single semi-circular notch in a plate was solved experimentally 
by CoKER and FILon [1]?), and M. M. Frocut [2]; theoretically by F. G. Maun- 
SELL [3], T. IsHrpasti [4], E. WEINEL [5], C. B. Linc [6], and H. Upocucut [7], 
using polar or bipolar co-ordinates, while the problem regarding the stresses 
when an infinite number of semi-circular notches are placed so closely as to 
exercise an influence on one another was theoretically solved by C. WEBER [8], 
and M. HEetTEnv1 and T. D. Liu [9] by methods enabling approximately correct 
solutions to be obtained. 

Applying the method MAuNSELL adopted, this paper gives a theoretical 
solution which strictly satisfies the various boundary conditions of semi- 
circular notches, and the numerically determined stresses in the neighbourhood 
of the notches. 


Method of Analysis 


Consider a semi-infinite plate containing an infinite number of semi-circular 
notches and receiving a uniform tension T. The radii of the notches and the 
distances between them are designated as a and b. Denoting the centre of an 
arbitrary notch of notches O and the centre at an arbitrary distance of 
kb(k=+1,+ 2,...) from 0, O,, we take the rectangular co-ordinates (x, ¥) 
and (%,, yx) designated as origins O, O,. Polar co-ordinates (r, 0) and (r,, 0;) 
are also used. They are assumed as shown in Figure 1. The following equations ~ 
are obtained: 

x=rsind, y=rcos6, (La) 
and 


Ke=x+kb, y=y, T= {(% + kd) + yh, tang, = =="? (1b) 


1) Department of Mechanical Engineering, Téhoku University. 
2) Numbers in brackets refer to References, page 477. 
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When the plate is in a condition of generalized plane stress, the aver age 
stress across its thickness can be derived from a stress function y which satisfies 


it Meh (2) 
If no notches are present, the uniform tension in the x-direction is given by 


r= 4 TP (1 + cos2 6). (3) 


MY 


See ae. 


Figure 1 


Notched plate under tension. 


To complete the solution when the notches are present, we first construct a 
series of functions which satisfies equation (2) and gives no stresses on the 
straight boundary and the same traction on each notch. With these functions 
added to y and all the boundary conditions at the rims of the notches satisfied, 
all requirements of the solution are fulfilled. This series of functions is derived 
by differentiating two fundamental functions which have Staeanitiee at the 
centre of each notch. 
As the first fundamental function, we take 


%o = —7 O sind — Y’r, 6, sin O, (R= +1,+2,...), (4) 
: 


which gives a normal force acting at the centre of each notch. 
To express this function in terms of the polar co-ordinates (7, 8), we intro- 
duce a complex co-ordinate @ referred to the origin O as pole 


Seah 26) 
From equations (1) and (5), we get 


log (9 + k 6) = logy, +74 (F — 0x) sa (6a) 
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and eis) ipa 8 
log (9 + k 6) = logh b+ 37 (325): (6b) 


From (6a) and (6b) the expansion of 6; is obtained as follows 
es 


AON Cos 2 Wak) oop Pye pa 


In the above expansion, some of the constant terms are omitted as they do 


not contribute to the stresses. 
By substituting equations (1) and (7) for equation (4), it is easy to deter- 


mine Yo: 


—)rrt 


. ( 
6, =» n+ 1) (k b)2”+1 


Xo = =r sind + 3 Sgn ranstleos(2m—1) 0+ 5 noi (2 +1) ot, 
(8) 
where the trivial terms are omitted and o,, is defined as follows 
Ov oS? init Tee eer ee (9) 
1 
Then yp gives stresses 
Rinne Thay ot way OS (5) vai 
oe pid ERMINE or ada ee DB pan O2n" 


n=1 


x {(2  — 3) cos(2m —1) 6+ (2 —1) cos(2n+1)6}, 


Oy (Han 
Cone bre a pen ORT hae 
a (10) 


x {(2 +1) cos(2 m —1) 6 + (2 —1) cos(2” +1) 6}, 
Od ro 7 (—)” pn 
T9= Bear (~ _ sf ie 02,7 *(2n —1) 


x {sin (2 n —1) 0 + sin (2m-+1)6}. j 


Therefore it indicates no tractions on the straight edge but is analytic in the 
plane except at the centre of the notch. 

Differentiation of y) with respect to x also gives functions with the required 
properties, but the odd derivatives must be excluded, since they are not even 
in x and cannot be used in the symmetrical system required. y, itself must also 
be excluded because it represents an unbalanced force which will have to be 


EE 
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balanced at infinity. The following series of functions is easily found, 


st setae B cya 
Kom (2 m — 2)! Ox2m y2m—1 
x {(2 m + 1) cos (2 m="1) 6 + (2 m —1) cos(2m+1) 0} 
= ee). 2 Le (2 — 1)! G2 2n—2m+1 alas 
eng (2m — 2)\(2n—2m+ 1)! 52" ” ™) 


x {(2 2 — 2m —1) cos(2m—2m+1) 60+ (2n—2m+1) 
x cos(2n — 2m —1) 6} (m = 1)’. 
As the second fundamental function, we take 


xo = —r 0 cosd = Dire Ge CosO ys, (12) 
: 


which gives a tangential force acting at the centre of each of the notches. It 
may be expanded in the same way as mentioned above, into the following form 


fo = —1 0cos0 + YC) Sen pants fin (2m +1) 6+ sin(2n—1)6}. (13) 
n=1 ~~ 


Here only the odd derivatives are required and the results are as follows 


- 2 (—)m+1 g2mr1 y’ 1 
X2m — 5 a inEE — “yam {cos2 m 9 + cos(2 m + 2) 6} 


< a os (PF 42 me =A)! O55, y2n—2m (14) 


(2m)! (2m — 2m — 1)! b?” 


n=m+1 
x {cos (2m — 2m) 0+ cos(2n—2m—2)0} (m 20): 


Now the required function is constructed in the form 
=< rar (1+ 0052.6) + 37 Anton + 37 Ba tit (5) 
m=1 m=0 


and then the problem reduced to determining A,, and B,, so that o, and 1, 
_ may be zero when s is a. 

Since o, and t,, contain terms of cosines and sines of odd multiples of 0 
respectively, it is necessary to expand these terms into Fourier series between 


6 = —x/2 and 0 = 2/2. The required expansions are 
lop es” (= mts (2m + 1) 
peosam+1)9=— amen +2 @msir—ast 828%, a 
s= | 1 


4 - 
sin(2m+1)0=— De ig koa sin2s 0. 
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Finally, we have the following expressions for the stresses 0, and T,9: 


[o.e) oO co 2 
( V2) (290 =i) a Ae 2m — On na 
phe ps; DRS pera 1(2m—2m-+1)!b™ 


m=1 n=m+1 s=0 

x (482? + 2m—2n—1)%"x, cos2 s 6 

bys ean m {m cos2 m 0 + (m + 2) cos(2m + 2) B} (17a) 
m=0 


Sy yy (=)" 2 (2m — 1)! a? 29? oon Bry 
iy MBP) Memes ich 


(2 — 2m) cos(2n— 2m) 0+ (2n—2m-— 4) 


x cos(2 1 — 2m — 2)6}4 = cos2 6) =0, 


“eA Soi ance ae 
ee ae (2m —2)!(2n—2m +1)! 62" 


Mesias. O 


2S cs dae ™ {msin2m 6+ (m+ 1) sin(2m-+ 2) 0} ¢ (17b) 


lo. <) co 
(—)"+12 (2m —1)! a2"-2m-25,, BL 
pan ON 2m)! (2n—2m—2)1 be 


+ 
(2 — 2m) sin(2n —2m) 6+ (2n—2m-— 2) 


x sin(2m — 2m — 2) 0} + > sin2m6=0, 


where |: 
my — ca be a 
0 2 (2m—1)? (2m +1)?’ 
fe dpe if (ar 18 
{(2m —1)?— 45%} {(2m 41)? —4 3° Sa} 


vie 16 (2m —1)2m(2m4+1) A, 


m 1 ) 
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pe soar ‘ (—)n+m off, 
6 2(2n—2m—1)*(2n—2m+4 1)?’ 
—\nr+m+s 
ty, = (—) 


{(2 — 2m — 1)? =4 52} {(2 0 — 2m +1)? — 4 53} ; 
ae 16 (2 — 2m — 1) (2m—2m) (2n—2m+1) Am 
m a : 


M13 
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(18b) 


The boundary conditions at the rims of the notches are thus satisfied if the 
coefficient of each term vanishes identically. After some reductions, we obtain: 


ss ae oy i=)" (2 Line 1)! eaten 6 Fon Arn (4 Wi 47 6) yes 


(2m —2m +1)! (2m — 2)! 52” 


+ ay (=)? 8 (2 x—2 ams2 Bm [ 
m=0 

J xied oe (2m + Sa canst Be 7 
m=0 

54S Ah m—1my— 2 

ia Ae een ae 

ie So Seed Si 3 wha? ea 
m=0 

+ yy et erin Be 0, | 
m=0 

& gith (2s —2m+1) a, — ee 


eis Lee) A 4 mk Fld 2) Pe, 


(2m — 2m-+1)! (2 m— 2)! b2” 


y(t" gs s (2 s+ 2m — 1)f a?*=? 0,2 am Bis 
Dy (2 ah (2 g.— 1) Gets 


(2 m)! (2 s)! h2s+2mr2 


» (srt 4 (2s +2 + 1)) a? de5 ame Bm 0 


(19a) 


(19b) 
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, 

co fe,2) (3) (2 ihe 1)! q2n-—2m—1 bya Ae - Ps 
Pe Ps, (2n—2m +1)! (2m — 2)! pin A oe oS ee . 
(19d 


m=1 n=m+1 


=e sy (—)*+™ AN Ques =i 2m ae 1)! Ca F2542m Bin j 
(2 m)! (2s — 2) beerem 


co 
(-)Mtt2 85. (29 46 2 m 4 1) lat? o x Boers 
) 2, ke (2 m)! (2 Shileimaaes = 0 (sie | 


B, can now be eliminated, and the following equations are obtained: 


1 Ba eo 


. ~ ai 2 (2 Rh 1)! qzn—2m—1 Oo ma At 
age a (2 —2m +1)! (2m— D1 ban (2m —2n—1)”™KXq 


p2m+2 =) 


—A 


1 Ae 
o> aaaer {(L — 29m) a, + G — 2m) a5} 


aes $y (—)" (2 wn — 1) at 2S Gg Ae 


m=1 n=m+1 2 n—2m aR 1)! (2 m — 2)! b2n 


“>> (2)! (2 ae 21) Cope Ban (20b) . 


b2m+2 


(—)™ 4 (4m + 3)! a oom Be 
i (@m)lbmre 


co 
(—)™4+1 (2m +5)! a4 Oomu, B 
a3 a 31(2m)| bamee m=—=0Q, 
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co A! 
D, gamit {(2 5 — 2m — 1) a+ (2s ~2 eX 3) ee 


m=1 


ey » (— ie (2 % —1)! ans is Can Ait 
(2 —2m +1)! (2m —2)1 62" 


m=1 n=m+1 


oe { (4 st4n—4m—2)"™a,+ (4s+4n—4m-+ 6) Pee Aci 


as Spc \m-+s as (2 m+2s'— At qrs-3 Oom +23 Bm (20c) 
TO (2 m)\ (2s — a 2)! p2m+2s d 

de hs (—)™t84! 8 kam + 2s +1)! a a2 ee B,, 
= (my! (2s) Bama By 

a (jets 4 (2m + 2s 43)Fa8 4? tices a Ben 

T 1 p2am+2s+4 m — (0) (s > 1) 
m=0 (2 m) | (2 Ss ae 2)! b 


A formal solution of the set in equations (19) and (20) can be found in Linc’s 
paper [10]: By assuming the terms including o,/b” as zero, the first approxi- 
mation is obtained. This is thus equivalent to MAUNSELL’s solution [3]; except 
that, in one pair of the corresponding parametric coefficients, the signs differ 
from each other. Then, with the coefficients determined in the first approxi- 
mation, the values of the terms including o,,/b” are calculated. As the constants 
of (19) and (20) this values being taken, the second approximation is deter- 
mined. Such approximations are repeated until the values given are convergent. 
Thus the values of A,, and B,, are obtained. Therefore, the approximations 
carried on subsequent to the first one are indeed corrections of the effects of 
the neighbouring notches. 

For brevity, no proof of convergence is pursued here, but it will be illustrated 
by the numerical examples which follow. 


Numerical Examples 


The foregoing solution was worked out for the cases b/a = 4, 6, and 10, with 
a as unity. The values of A,, and B,, determined are shown in Table 1. On 
working out the values of o, and t,, at the rim of the notch, by using the 
values shown in Table 1, the greatest values found were 0-003 7, 0-013 T, and 
0-018 T for b/a = 4,6, and 10 respectively. These are the values of o, for 
6 =2/2. Therefore the values shown in Table 1 are thought to be sufficient to 
get the stresses with a high degree of approximation. 

The values of 
o, on the edge of the semi-circular notch (Table 2), 
o, on the straight edge (Table 3), 
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a 


9—O1T X 12 9-OTX. o-T— 
e—OTX L9:2 9— OT X 9-T e—OLX T6-T 9— OT X cop 
eOLX 86-7— 7! 9—O0T X SL | e— ODA S656 Sao Oles Reb 
e—OTX T0-8 9—OT Xx 0-2 e—OTX 62:9 9—OT X oo 
p—OLX 8Z€-T— | 9—OTX LBs =| sy = 0F 8505 Feo 07 VE 
y—OTX S022 9—OT X Cl Wee abes Yara 9—OT Xx or 
y—-OTX vOL-€— | e—OTX 82-T y—OTX CET-E— | 9—OTX LG 
+—OT X 888-9 e—-OTX 96:S5— | s-OTX Z9L-¢ oe Ou) re 
e—OLX pSE-T— | s-OTX 6SL-2 e—OLX O9T-T— | s-OTX O8L:T 
e—OTX OZO-€ e-OTX L9:-T— | e—OLX [06-2 OES CYA Ea]t 
e—OTX TS8-8— | z—OT X T+H09-2 z—-OLXLPLT-T— | g—OL X009b-2 
1-OE XLL86-T — | 1—OT X9800-€ 1-OT XST08-T — | 1-OT XTET9-Z 
t— OT XZ88T-2L — TOT X6OT12:9\— 

aie We Ub ye Wg | Wy, 
Or = »/9 = »/9 
fe On 


g “py squats fao9 
T PGP 


9—OTX sb ot 
9—-OT* Ly 9-OI Xx Le IT 
g-OTX $8:T— | 9-OTX 9-0 — Or 
e-OTX 86-2 9-OT* ov 6 
o—OTX FeS— |} 9-OLTX 6:9— 8 
s-OLX PST-T e-OTX SOT L 
s—OTX SST-C— | g-OTX 66-1 — 9 
v-OTX E7-S g-OLX E72 ¢ 
s-OTX ¥OT-T—-| e-OTX O€4%— v 
e—-OLX €86-€ y—-OTX SLL-4— € 
z—OEX€crO-T— | z-OLXTELS-2 c 
1-OT X 8S69-:T— | 1-OT X PLPT-2 ib 
1-OT X SOTO-S — 0 


y= »v/9 
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Table 2 Table 3 


o,/I for r=a o,/T for 0 = 90° 

6° | ba =4 | va=6 | b/a=10 +16 

0 2-13 2-46 2:79 1-0 
15 1-91 2-25 2:56 11 
30 1:34 1-67 1:94 1.25 
45 0-65 0-94 1-13 155 
60 0-12 0:30 | 0:39 1:75 
75 —0:06 | —002 | —0-01 2-0 
90 0-00 0-00 0-00 


Table 4 Table 5 
a, for = 0° Oy for 0 = 0° 
| rla_ | bla=4 | ba=6 | dja = 10 rla | b/a=4 bla=6 | bla = 10 


| 1-0 0 0 0 1:0 2:13, + 


2:46 2°79 
Hig O15 0-18 0-20 ike ds f3 2:01 Zs 
E25 0-23 0-29 0-33 1525 1-40 1-61 1-82 
| Ne 0-24 0-31 0:37 | io) 1-14 1-28 1-43 
1-75 0-20 0-28 0-35 ihe 1-04 Luc) 1-25 
2°0 0-16 O24" =| — 0-30 2-0 1-00 1-05 1-14 


o, and og along the symmetrical line Oy (Tables 4, 5), have been calculated, 
and are here given. 

The results of Tables 2-5 for b/a = 4 are shown in Figure 2, where the 
corresponding stresses given by MAUNSELL are also plotted for comparison. 

The striking features of the results are that the maximum concentration 
of stress becomes smaller, and the tension on the straight edge changes into 
compression, though small in amount, as the neighbouring notches are placed 
more closely together. The reduction of the stress concentration is shown in 
Figure 3. 

It will be noted that the maximum concentration of stress for the case of 
b/a = 4 is 2:13 T, which is just coincident with the result given by C. WEBER[8], 
while that given by M. HetEny1 [9] is 2-175 T. 
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Figure 2 
Stresses in the neighbourhood of each of the notches. 


I stresses at the rim (Table 2); ZZ radial stress on the straight edge (Table 3); III radial stress 
on the symmetrical line (Table 4); JV tension across the symmetrical line (Table 5). 
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Zusammenfassung 


Es wird hier ein Lésungsverfahren fiir das Spannungsproblem der Halbebene 
periodischen Halbkreiskerben nach Figur 1 besprochen. Zur Erfillung der 


Randbedingungen werden Spannungsfunktionen mit Singularitaten, die perio- 
disch zur x- und symmetrisch zur y-Achse sind, in den Mittelpunkten der Halb- 
kreiskerben angenommen. 


der 


Die Spannungserhéhung im Kerbgrunde und auch der Spannungszustand in 
Nahe des Randes werden zahlenmassig bestimmt. 


(Received: February 16, 1957.) 
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On a Problem of Heat Conduction with Time-Dependent 
Boundary Conditions 
By TsE-Sun Cuow, Warren, Michigan, USA?) 


1. Introduction 


In a recent paper [1]2), Vopi¢éKa gives general results for heat waves in a 
homogeneous, isotropic cylinder whose boundaries are subjected to periodic 
temperature fluctuations. The purpose of this paper is to derive the results for 
a more general and complete’problem: the determination of the distribution of 
temperature in a homogeneous, isotropic cylinder when its initial temperature 
distribution is known and its boundaries are subjected to temperature fields 
which are arbitrary functions of time. The method we shall use consists of a 
substitution which reduces the non-homogeneous time-dependent boundary — 
conditions to homogeneous boundary conditions, and may be considered as 
an extension of the method used by MINDLIN and GooDMAN in solving vibra- 
tion problems with time-dependent boundary conditions [2]. The method is, — 
however, quite general, and may be applied to problems of similar types. 


2. Statement of the Problem and Removal of Non-Homogeneous 
Boundary Conditions 


Consider a homogeneous isotropic cylinder, 0 < @ S7, |z| Ss subjected 
to time-dependent boundary conditions and an initial temperature distribution 
at t= 0. It is desired to find the temperature distribution in the cylinder at 
some ¢ > 0. Expressed otherwise, we are required to find u(@, , z; t) from the 
following equation: 


Ou ae 1 ou 1 072u 0?u | 

past a fel : 2 : 

ot a 

| 00 @ 00. 03 0G dea (2.1) 
(Oss O'S %, OS Si Aare Vices 


with the following boundary conditions 


et hy u= fp, 2) a(t) (o=7), ot + hh w= fe (0, ) a(t) (z=s), 


ae (2.2) 
Cryo wer fs (0, p) ast) (z= —s), 
1) Research Staff, General Motors Corporation. — The author wishes to thank the Cor- 


poration for permission to publish this paper. 
2) Numbers in brackets refer to References, page 484. 
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where fy, ha, hs > 0 and the initial condition: 


u(o, P, 2; t) = ul, 9, 2) 


By putting 


u(o, @, 2; t) 


f= 0). 


3 


l=1 


= plop, 2:1) + J” Xi(o, @, 2) Ti(t) 
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(2.3) 


(2.4) 


and substituting this expression into (2.1), (2.2), (2.3) we find the equations 


for p(o, y, z; t) are 


oy 2 { O° 1 oy i 

Bee at Zant 
3 

02x Ox 1 

= q2 eel] Z LOANS peed 


0 

5 +h y= Aly, 2) alt 

oe thay = ble, 9) 
- 2) = [alQ, P} 


Ls hs yp = f5(0, p) &(t) — 


WO, Y, 23 


We now choose 7;(¢) as follows: 


es eae 
0g? Oz? | 

2X OX py 
0g? ez2{.-F 


hs x} T, (z 


0) = u% — 


Ns X,T,(0) . 


l=1 


ee ai (2 


I 


T(t) = a(¢) , 
while X, are to be determined by 
ox, 1 Ox, | 1 etiXe 08X, 1 
do? + 9° 09 1 oF Og? * dz? 


with the boundary conditions: 


As + hy X,= 61, (9, 2) (o=7), 
ae + he X_= Oyo f(0,p) (2=5), 
GED hg X,= brs fle. 9) (2 = -8) 


where 6, = 1 for / = m and is otherwise zero. 


ea dys 
1 


| es 
| 


(2.6) 


(2.8) 


(2.9) 


(2.10) 
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With this choice of X, and T, we find the equations that determine p are 
given by: 


oy 0?y Lie Ow J ial had Ory | XT (2.11 
i. oe ae + 0 00 ie ie ee “Oz? | --¥ U Ty, ) 
) 0 
St + hy =0 C=); E+ hey =0 (z=s), 
. (2.12) 
oy 
aE Tig Oat ale 
3 
v(o, p, 2; 0) = % — J) X, T,(0) - (2.13) 


Thus, when X, and T;, are properly chosen, the problem of determining u is — 


reduced to one of determining X, satisfying (2.9) and (2.10) and yw satisfying 


(2.11) with the homogeneous boundary conditions (2. 12) and the initial condition 
(2.13). We proceed first to find X;. 


3. Determination of X,(0, P, Z) 


We notice that each X, satisfies the differential equation (2.9) and two } 
homogeneous boundary conditions and a non-homogeneous one (2.10). Follow- 
ing the classical method of separation of variables we expand each X; into a — 


doubly infinite sum of the form R(e) P(g) Z(z). We present the final results 


here, omitting the necessary details. It is suggested that readers refer to [1] © 


for the proof of the orthogonality property of the characteristic functions and 
the evaluation of the normalization integrals. 


co a 


(0, P 2 me yids o)e en? Loos S24; sin 2} G sve 
ak (3.1) 
where 4,;, €;, C; and ,A,,; are given by 
(hy hs s” — 6?) sin2 6; + (hy + hg) s 6; cos2 6; = 0, 


, 0; 7 


€, = —_, C= 


6; sin 6; — he s Cos 6; 
0; CoS 0; + hy s sind; ’ 


+s 
, &; : ’ ei ai . 
An j \( ip hy) [,,(é;) ae = Ih 4x(€) | {cos z+C; sin 2 z\ dz (3.2) 
+7 


-ae/ eos % 2+ ¢, sin % aja | f(y, 2) e~'"? dp. 
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= = € On k oO; 
X,(0, y, 2) = py ay; Ee f, (“28 g) ¢ Pe {cos -2-+ Ci, sin—* E al (3.3) 
k=1 n=—co 
where On%; Eng Cry and 24, are given by 
(n + hy r) Tn (En x) + Enk Tn+1(Enx) =0, 


6’, — fnks Cc! hs S COS Onn — On nk SIN Ong 
nk AS nk 


~~ , 


¥ 


5008 645 Ay S SIND « ” 


a Onk S (Re + hs) Cos2 Ong + (S* hg hg — Onn) SINZ On (3.4) 
Set $ (dg; cos 0, + Ay 5 sin d;,) 


x [of i( o)} dg ~ai fet “(2 ode [ He g) ei"? dp, 


X5(0, p, 2) -> “> Ane I,(“2 


=1 =—©o 


) ef” Joos Pat z+C),sin- One ae ah, (3.5) 


where i, w Eng Cnz and 3A, are given by 


(nm Bo hy r) T(En x) + Enk Tn +1(Enx) =0, 


r _ EnkS cr Ong SIN Of, — Ng S COSO/ 
wh ieee nh ~ “61, COS 64, + hg $SiN Oy,” 
34 * Ong S s (hg + h 3) cos 2 Onk + (s? hs hs -e On x) sin 2 Onk. : (3.6) 
n 


ida Gos 025 A, © ein 02) 


Yr +76 
nk 2 1 n —tn 
x fo{z(2* e)| do=s— fot (=2* g ) [tle 9) ¢ dg. 
0 0 —71 


4. Determination of y and u 


The function yp satisfies the differential equation (2.11), the homogeneous 
boundary conditions (2.12) and the initial condition (2.13). This suggests that 
wy can be expanded into the triply infinite sum: 


yo, p25 4) = 35) S F, x s(é) I," o) 2"? {cos z+ €, sin zh. (4.1) 


From the boundary conditions (2.12) we see that ¢,,, 0;, C; are determined 
as follows: 


(n + hy 7) In(Enz) + Ent IngalEnz) = 9, | 
(hy hg s* — 6?) sin2 6; + (he + hs) s 6; cos2 6; = 0, (4.2) 


6; sin 6; — hy s cos 0; 
6; cos d; + hy s sind; 


¢;= 


_ ZAMP VIII/31 
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It remains only to determine the coefficients in (4.1), ie., Frx,(/). To this 
end we substitute (4.1) into (2.11). This results 


where use has been made of the fact that J,(€,; 0/7) satisfies the following 
equation: 


at ieee 0) ae . . i 1, ote 0) (22 oa 7) Tal ae 0) =0. (4.4) 


Wey y? 0 Y 


To solve for F,;,,(¢) in (4.3) we expand X,(o, y, 2) into the following triply 
infinite sums: 


Xie eke Eh By Ee 1,( se 0) eit? {cost 2 + C; sin a}. (4.5) | 


Using the orthogonality property of I,(én,%0/r) and {cos (0,;/s) z+ C, sin (6,/s) Pe 


we find the relations between ,B,;,; and ,A,,;,; to be, from (3.1), (3.3) and (3.5): — 


asf etl ye) de stuns 4( Se) ol 0) do, 


+5 
B Py ara 
2Bnng | \COS-™ z+ C; sin-* 2 dz 
+s 


Onn 
= oAnki {cosa = Pz+C),si 


=—s 


+s 
Se ie Se 
Bars | {cos 2 + Cicin z} dz 
—s 
+s 


5 
= Anns {cos = 24 Cl sm oe z} {cos z+ C,sin 24 z} dz. 
Ss Ss 


—S) 


By substituting (4.5) into (4.3) we obtain by ua coefficients : 


, Enk 6% 
Fin!) eke en oo, | i mki(t) = -3 1Ba ng % (¢ (4.7) 


al {cos z+ C; sin °2 z| dz, (4.6) — 
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and upon integration, we have immediately 


2 2 
enk oF 
GO} et piace Fg 
y? oh 
e 


Fing(t) = Fax3(0) exp 


3 “8 ; € 82 
nk ; n ; 
ws Bn kj exp] a® = — =~) i} | a, (T) exp} ( is — S r| at, 
0 
1. G.; 
2 nh i \ 
Fx j(t) F, ¢3(0) exp{ a ( Ty ie 4 a | 


saan fuk SF tie 
$0 (% A) avonlo (4-24 


- f é? 
a, (t) exp| a2(- nk a 
U - ye s2 

0 


by integration by parts. The coefficients F,;,;(0) are determined by the initial 
condition (2.7). Let u9(0, y, z) be expanded into the triply infinite sum: 


Uo(O, P, 2 ee Ds Dye; 1,(“2* @) 2 Joos 2 + C;sin 2h, (4.9) 


then from (2.7), (4.1), (4.5) and (4.9) we obtain 


3 
Fix (0) i" Dy xj =) 08 nk 9g a,(0 (4.10) 
F=0: 


From (4.1), (4.8) and (4.10) we therefore obtain 


— 


v0, P, ; t) 


3 ee, 62 Enh Of 
—a S- i e +) Bras €XP|? (3 e 3) i (4.11) 
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3 
Remembering that u(o, y, 2; 1?) = pla, @ 2; 4) rT ges X,(0, Y, 2) a(t) and using 
(4.5) we finally obtain from (4.11) the solution in the following form: 


(9, 9, 2; t) | 
& opp S) Dans exp|a? ee a a | 
A ay (ee _ -) 1Byxs exp{a? (3 al a) | (4.12) 
e fete exp| a? (ss 2 *) t| me 
: . 


cae (“2 0) ZEe. {cos.% z+C; sin zh, 


if 


where ,B,x;, Dax; are given by (4.6), (4.9) respectively, and ¢,,, 0;, C; are given 
by (4.2). 
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Zusammenfassung 


Das folgende Warmeleitungsproblem wird analytisch behandelt. Auf dem 
Mantel und den Deckflachen eines Kreiszylinders sei eine zeitlich verdnderliche 
Temperaturverteilung und ausserdem die Temperatur einer Kurve zur Zeit t = 0 
gegeben. Gesucht ist die Innentemperatur fiir alle Zeiten. Die verwendete 
Separationsmethode kann auch auf andere Randwertprobleme angewendet 
werden. 
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Design of Circular Plates for Minimum Weight!) 


Emin TuRAN ONaT?), WALTER SCHUMANN®) and RicH. THORPE SHIELD?), 
Providence, R.I., USA 


1. Introduction 


The minimum weight design of circular plates has been investigated by 
Hopkins and PRAGER [1]*) and FREIBERGER and TEKINALP [2] for the case of 
a uniformly loaded plate with a simply supported edge. Intuitive arguments 
were used in [1] to obtain the minimum weight design, and [2] employed the 
calculus of variations. In this paper the direct design procedures developed by 
DRUCKER and SHIELD [3, 4] are used. These formal design procedures were 
derived from the theorems of limit analysis [6]. Both built-in and simply sup- 
ported edge conditions are considered and the plate is loaded by an arbitrary 
rotationally symmetric distribution of pressure. The designs obtained are abso- 
lute minimum weight designs for the sandwich plate but the minimum is a 
relative one for the solid or homogeneous plate. 


2. Statement of the Problem 


The problem under consideration is the design of a circular plate under 
rotationally symmetric types of loading and edge support. The design is to be 
such that the plate is just at the point of collapse under the given loads and such 
that minimum volume of a given elastic-perfectly plastic material is involved. 
The material is assumed to be homogeneous so that the design for minimum 
volume coincides with the design for minimum weight. As in work on the load- 
carrying capacities of circular plates (see [5], for example), the term collapse 
is used in the sense of limit analysis [6], and membrane stresses are ignored in 
the analysis. 

The plate is assumed to be horizontal and loaded by a distribution of 
pressure (7) on its upper surface, y measuring distance from the plate center. 
The edge ry = R of the plate is simply supported or built-in. The radial and 
circumferential bending moments are denoted by M and N, respectively, with 
the convention that positive values of these moments stress the lower surface 


1) The results presented in this paper were obtained in the course of research work sponsored 
by Watertown Arsenal Laboratory under Contract DA-19-020-ORD-2598-DA-19-020-501-ORD- 
3984. 

2) Brown University. 

3) ETH and Brown University. 

4) Numbers in brackets refer to References, page 499. 
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of the plate in tension. For equilibrium the moments satisfy the equation 


d\ryM) aN i. 
oN trp =o. (1) 


The shear force Q is given by 
rQ=—,—_ _- NN. (2) 


The curvature rates x, Ain the radial and circumferential directions are given by 


A d 
(a aOR 3 
where w is the deflection rate of the middle surface of the plate in the down- 
wards direction. 

For simplicity it is assumed that the elastic-perfectly plastic material obeys 
TreEsca’s yield condition and the associated flow rule. In terms of the moments 
M and N, plastic yielding can occur when one (or both) of the magnitudes of 
M and N equals the fully plastic moment M, or when the magnitude of M— N 
equals M,. The condition defines a yield curve in a plane in which the moments 
M, N are taken to be rectangular cartesian co-ordinates and the yield curve is a 
hexagon (Figure 1). For a point on a side of the yield hexagon, the flow rule 


Figure 1 
Moment yield condition. 


associated with the yield condition requires that during collapse the curvature- 
rate vector (x, A) is parallel to the outwards-drawn normal to the sides. At the 
corner of the hexagon, the vector must lie between the normals to the side 
which meet at the corner. 


The rate of dissipation of energy D, due to plastic action per unit area of 
the middle surface of the plate is given by 


De=Mxnx+Na. (4) 
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Because of the flow rule, the rate of dissipation D, is uniquely determined by 
the curvature rates x, A. 

In the problems considered here, replacing the Tresca yield condition by 
the perhaps more acceptable condition of von MisEs would make an analytic 
solution impractical and a numerical approach would prove necessary. A simi- 
lar situation exists in the previously mentioned work on the load-carrying capa- 
cities of circular plates. However, when an approximate yield condition is used 
to estimate load-carrying capacities, bounds to the error committed can be 
determined [7]. As will be shown, an analogous procedure applies to minimum 
weight design. For the purposes of illustration we shall assume the von Mises 
yield condition with yield stress o) in tension to be the true yield condition for 
the material and the Tresca yield condition to be the approximate condition 
used to simplify the problem. The extension to any true yield condition and 
any approximation is obvious. The yield surface defined in principal stress 
space by the von Mises condition with yield stress og circumscribes the yield 
surface defined in the same space by the Tresca condition with yield stress dp. 
It follows that any structure designed to be just at the point of collapse under 
given loads for the inscribed Tresca condition will be safe or at collapse under 
the same loads for the von Mises condition. Thus the volume V7(o,) of the 
minimum volume design based on the inscribed Tresca condition with yield 
stress oy will not be less than the volume Vy(o9) of the minimum volume 
design based on the von Mises condition. By the same reasoning, the volume 
Vy(o,) of the minimum volume design based on the von Mises condition will 


not be less than the volume V7(2 6o/V3) of the minimum volume design based 


on the circumscribed Tresca condition with yield stress 2 0, /V3. Thus the true 
minimum volume Vj;(¢9) is bounded from below and above, 


V4 ve) S Vy (Fo) S Vr (oo) - (5) 

The accuracy to which Vj,(o,) is determined by (5) depends on the variation 
of V;(o,) with oy. In the sandwich plate examples considered in section 3 
below, V;(o,) is inversely proportional to og and (5) determines Vy(o9) to 
within 7%. In the solid or homogeneous plate examples of section 4 below, the 
designs obtained are relative minimum volume designs with volume inversely 
proportional to the square root of o». It seems probable, however, that the 
designs are in fact absolute minimum volume designs, and with this conjecture 
(5) determines Vy,(aq) to within 4%. 


3. The Sandwich Plate 


The sandwich plate has a light-weight core of constant thickness H between 
two identical face sheets of variable thickness ¢ («<< H) composed of the given 
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material. The minimum volume design is just at collapse under the pressure 
p(r) and involves minimum volume of the face sheets. The core carries no 
bending stresses and the fully plastic moment M, has the value 


Mi iy tie (6) 


where oy is the yield stress of the cover plate material. 

It has been shown [2, 3, 4] that a sandwich plate designed to collapse in a 
mode such that the rate of dissipation per unit volume in the faces of the plate 
is constant is an absolute minimum volume design. In terms of the rate of 
dissipation D, per unit area of the middle surface, D,/¢ is constant for such 


a design, 
Dy Mxzxt+twNA 


t t 


= const . (7) 


This condition is independent of the thickness ¢ because of the linear depend- 
ence (6) of M4, and therefore of D, on ¢. With the Tresca yield condition 
(Figure 1), the condition (7) restricts the position of the moment point on the 
yield hexagon. For example the moment point cannot lie on the side AF for a 
finite range of 7 because the flow rule requires A = 0 and the condition (7) 
cannot be satisfied in view of (3). It is readily found that the condition (7) can 
only be satisfied for a finite range of y when the moment state is represented 
by one of the points A or D, C or F of the hexagon. The points A and C will be 
used in the following. 
For the point A, condition (7) becomes 


d?w 1 dw 


dr? se yy ee ee (8) 

where « is a positive constant. It follows that 
ic: Bee Seema eae 9 
4 of I 2 a Obs ( ) 


where } and f are constants. In addition x and 4 arenon-negative at the point A 
and this enforces 7 = b. 
For the point C, the condition becomes 


d?w 
Va - (10) 
so that 
1 
Dis sot (2 Oa (11) 


where c and y are constants. The requirements — x =A =O restrict the 
range of y toc =r = c/2. 
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When the plate is simply supported at the edge y= R and loaded by a 
downwards-acting pressure (7), the assumption that the curvature rates x, A 
are both non-negative leads to the minimum volume design. The moment state 
through the plate is then represented by the point A in Figure 1. At the edge 
y = R the deflection rate w is zero and at the center dw/dr is zero since A is 
finite there. The appropriate form of (9) is therefore 

1 


Wm ym (R? — 7?) . (12) 


For the point 4, M = N = M, and substitution in the equilibrium equation (1) 


yields 
$ ra? M aM 
Sd at + = =—7 pir) : (13) 


At the edge of the plate M = 0 and the shearing force Q given by (2) is zero 
at y = 0, assuming that there is no concentrated load at the center. With these 
conditions, integration of (13) gives 
R § 
r 1 pr 
M, = 0 Ht =| = d&/ 0 ple) de, (14) 
0 


ra 


Y 


and the thickness ¢ given by (14) is the minimum volume design [since f(r) is 

non-negative, ¢ is non-negative}. A ring of uniform force per unit length on the 

circle y= a can be included by suitable interpretation of the integral with 

respect to 9 in (14). The particular case of a concentrated force P at the center 

adds a term pu te, R is) 
2% Oe 


to the right-hand side of (14), producing an infinite value of ¢ at the center but 
the face sheets have finite volume. The result for uniform pressure has been 
obtained previously [1]. Further examples of simple types of loading are 
studied in more detail below. 

For the plate with a built-in edge and non-negative f(r), near the edge of 
the plate where x <0, A = 0, the point C of Figure 1 will apply. Near the 
center of the plate x and / are non-negative and the point A applies. Assuming 
that regime A applies for 0 < 7 < aand regime C fora <7 < R, equations (9) 
and (11) govern the deflection rate in these regions respectively. At the built-in 
edge w and dw/dr are zero and at the center dw/dr is zero. Also w and dw/dr 
are continuous at the junction 7 = a of the two regions®). In order to satisfy the 


5) Tt will be seen that the thickness ¢ of the face sheets in the minimum volume design is zero 
at the junction y = a and collapse modes exist for which dw/dr is not continuous across 7 = a. As 
we require a collapse mode in which the rate of dissipation per unit volume is constant in order to 
show that the design is a minimum volume design, discontinuities in dw/dr are not permissible in 
this type of collapse mode. 
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conditions it is found that we must have a = 2 R/3, and we obtain finally 
ul 
ae 


[Ferm GG Rsrse), 


a (2 R® — 37?) (osrszR), 


v= (16) ‘ 
3 


For 0 <7 <2R/3, M=N=M, and M, satisfies equation (13). In the outer 
annulus 2 R/3 <7 < R, M = —M, and N = 0 so that for equilibrium 


sx (yr My) = 1007) - (17) 


At ry = 2 R/3, M is continuous and therefore zero and the shearing force Q is © 


continuous. As before Q is zero at 7 = 0 in the absence of a concentrated load 
at the center. With these conditions, equations (13) and (17) give on integration 


2R/3 f 
[aos (0o<rs7R), 

Mi =o,Hi= (18) 
1 2 
i at [0 10) de (;RsrsR) 
2R/3 0 


Rings of force per unit length can be included as for the simply supported plate. 
The case of a concentrated load at the center and other simple types of loading 
are examined below. 

The case when part of the plate is loaded by pressure on the lower surface, 
that is when (7) is negative for a range of 7, will not be considered here. In 
addition to the points A and C, moment states represented by the points F 


and D of the hexagon may have to be used, depending on the pressure distri- — 


bution p(7). 
(i) Circular Loading, Simply Supported Edge 


In this case £(7) has the constant value p for 0 <r <7 and is zero for 
% <7 ZR. Direct substitution in equation (14) gives the thickness ¢ of the 
minimum volume design, 


A il 
GPR 1) + 5 pr log> V=a7a7e 

OHt=) . : (19) 
OPI ae 0) =r een 


The volume V of the face sheets is found to be 


Ah r\ apr R? 
@ x (2 as aT (20) 
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As might be expected, the volume V for a given total load a # 7? decreases as 
the area over which the load is spread increases. It can be seen from (20) that 
the average thickness depends linearly on the total load a # 7? and otherwise 
only on the ratio 7)/R, for given o, and H. 

The thickness ¢, for a plate with face sheets of constant thickness which is 
just at the point of collapse under the same loading and support conditions can 
be deduced from the results of Hopxins and PRAGER [5]. The volume V, of the 
plate with constant face thickness is found to be 

4 2 2 
Aes as 
Comparison of (20) and (21) shows that the saving effected by the minimum 
volume design over the constant thickness design increases from 25% to 50% 
as Y)/R decreases from unity to zero. 

The case of a concentrated force P at the center of the plate can be obtained 
as a limiting case by letting a p 77 tend to P as 7, tends to zero. The thickness 
distribution is given by 


q,Ht= ns log : (22) 
and the volume of the face sheets is 
ik PR? 
V = icy . So H Fy (23) 


(ii) Circular Loading, Built-in Edge 


The form of the thickness distribution depends on whether the radius 7, of 
the loaded area is less than or greater than 2 R/3, the radius of the junction 
between the regimes A and C. 

For 7) < 2 R/3, equation (18) gives 


TPR + 5 ples (SrSr), | 
oHt= = blog (S71 <= R), (24) 
1 oR ¥— 2 R/3 (2 R<rsR) | 
For 7) > 2 R/3, the corresponding result is 
i p(S Rn) (0<7s= Rj, 
ee Ear, (25) 


3727 — 8 R3/27 — 24% 
Y 


1 
6 
1 
6 


p 
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The concentrated load P at the center is a particular case of (24) and the 
thickness is given by 


/ 2 
i log - (0 = x, R), 
Oo, Hti=: , (26) 
Bas Nearer tne ({R<rsR). | 
2G Y 3 7 han 
The volume V of the face sheets has the value 
1 Pie 
V= cay (27) 


As has been shown by Hopxrns and PRAGER, the limit value of a concentrated 


central load P is 2 M, for a circular plate with constant fully plastic bending ~ 


moment M, for both the simply supported and the built-in plates. Thus the 
face sheets of the constant thickness plate which is just at collapse under the 
central load P have the volume 

VeNGe 
Ve 6,H , 


(28) 


irrespective of the support conditions. The saving effected by the minimum — 


volume design over the constant thickness design is 67%, for the built-in plate | 


compared with 50° for the simply supported plate. 


. 
4 


The thickness distribution for the built-in, uniformly loaded plate follows . 


from (25) by giving 7, the value R. The minimum weight design for this case is 


shown diagrammatically in Figure 2(a). The volume of the face sheets has the ~ 


value 
ap R* 


V = Ole (29) q 


0-664 R ——+| 


Figure 2 
Minimum weight designs for built-in circular plate under uniform pressure: 
(a) sandwich plate, (b) solid plate. 


Ew a. eee 


Vol. VIII, 1957 Design of Circular Plates for Minimum Weight 493 


With the results of [5], the volume V, of the facesheets of the constant thickness 


plate is found to be 
ap R* 


V, = O.078 SP 


(30) 
so that the minimum volume design effects a saving of 34°, over the constant 
thickness plate. 


(i) Annular Loading 


The results for annular loading, that is #(r) has the constant value # in the’ 
annulus 7) <7 < +, and zero elsewhere, can be obtained from the general 
formulae (14) and (18). For both the simply supported and built-in plates, the 
thickness ¢ is constant in the central region 0 <7 <7. In particular, the 
thickness ¢ is zero in the unloaded central region for the built-in plate when 7, 
is greater than 2 R/3. 


4. The Solid Plate 


For the solid or homogeneous plate, the fully plastic moment M, is given by 
1 
MM, = “4 % h*, (31) 


where / is the (variable) thickness of the plate. The results are less definite and 
the analysis is less tractable for the solid plate than for the sandwich plate. As 
has been shown [2, 3, 4], a solid plate designed to collapse in a mode such that 
the rate of dissipation per unit volume is constant on the surfaces of the plate 
is a relative minimum volume design. Here ‘relative minimum’ is used in the 
following sense. If h,, is the thickness distribution of the relative minimum 
volume design, then any neighboring plate with thickness h = h,, + 6h, where 
| dh | <<h,,, which is safe or just at collapse under the given loads will involve 
a greater volume of material. In terms of the rate of dissipation D, per unit 
area of the middle surface, the condition for a relative minimum volume 
design is D, ee a 


Se teas ; = COnSis (32) 


Since D , is proportional to h?, this condition involves the design thickness h, a 
circumstance which complicates the analysis. Thus in contrast to the sandwich 
plate, it is not immediately apparent that points A and D, C and F of the yield 
hexagon in Figure 1 are the only moment states involved in a relative minimum 
design. In the following, moment states represented by the points A and C 
will be used to obtain solutions. 

Considering first the simply supported circular plate, loaded by a non- 
negative pressure #(7), the use of moment states represented by the point 4 
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provide the relative minimum volume design. The assumption M = N= M, 
and the equilibrium equation together with the conditions at the center and 
edge of the plate lead to 


R § 
My = _ ooh? = |= dE | 0 Plo) de. (33) 
0 


£ 


as in the previous section. For non-negative p(7), 4” and therefore / is a mono- 
tonic decreasing function of 7. In order to show that the design thickness h 
given by (33) is a relative minimum volume design, it is necessary to show that 
an associated collapse mode exists satisfying (32). For the point A, (32) becomes 


d?w 1 dw ea 


dr? e Bae iy, | clei das (34) 
where « is a positive constant. As dw/dr is zero at r = 0, it follows that 
dw 4 0 
Ghee: -$ lay ae Go 
0 
and therefore since w is zero at the edger = R, 
R é 
= 1 Q 
w af = dé [=e do . (36) 
r 0 


In order that the collapse mode be associated with the point A, dw/dr and 
d?w/dr® must be non-positive. It follows immediately from (35), 


and since h(r) is monotonic decreasing 


d?w ea ch a 


SHG | O88 Fey = — FFG 39. (37) 


Hence the thickness h given by (33) provides a relative minimum design. 
For the built-in plate under pressure £(r), we assume as in the previous 
section that regime A applies for0 <r <a and regime C fora <+< R.In 
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the same way that equation (18) was derived we obtain 


[ra fopde (<r: 


0 


= 
| 
S 
I 
WA 
SS 


(38) 


r g 


1 fas fo ple) de (asr 


a 0 


IA 
& 


The thickness h is zero at y = a and for non-negative f(r) is monotonic decreasing 
for 0 = ry < aand is monotonic increasing fora <r < R. It remains to deter- 
mine a and the associated collapse mode w satisfying (32). ForO <r <a, w 
satisfies equation (34), where h is now given by (38), and as before it can be 
shown that dw/dr and d?w/dr? are non-positive in this region. Fora<7< R, 


condition (32) requires 
d?w od 


ees 39) 


Since dw/dr is zero at the built-in edge it follows that 
d ‘ d 
jw + Q 
= a“ fa (a<r<R). (40) 
At the point C, —x =A = 0 and with (3), (39) and (40) this requires 
BP BSD: POE Rey Sy eR (41) 


The radius a is determined by the condition that dw/dr is continuous at 7 = a. 
For 0 <7 <4, dw/dris given by (35) and with (40) the condition becomes 


a R 
ee. 3 ah ees 
J hoe a ie). ia) 


_ When a has been determined from this equation for a given p(7), the inequality 
(41) must be checked. Inequality (41) is always satisfied if a 2 R/2. In the 
range a <7 <R, h(r) is a monotonic increasing function of 7 and it follows 


that in this range 4 


rd We = - 1) a ; (43) 


Y 


Comparison of (41) and (43) shows that (41) is satisfied if a = R/2. 
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(i) Circular Loading, Simply Supported Edge 


When the plate is loaded by a constant pressure # over the area0 S7 S7%, 
the thickness distribution for a simply supported plate is obtained by direct 
substitution in (33), 


D) 


p \3 | ea 
(Ly fa + 2 


(44) 


(Enter ses 


The result for uniform loading, 7) = R, has been obtained previously by Hop- 
KINS and PRAGER [1]. The volume of the plate is proportional to the square 
root of the total load x # 7? and inversely proportional to the square root of do. 
As the ratio 7,/R decreases from unity to zero the saving in volume effected by 
the design (44) over the constant thickness design increases from 18° to 37%. 

For a plate with a concentrated load P at the center, the relative minimum 
volume design is 


jee (z-)" ( 2 log =) (45) 
with volume 
Va5 (=)" mR. (46) 


(i) Circular Loading, Built-in Edge 
As for the sandwich plate, the form of the thickness distribution depends 


on whether the radius 7, of the loaded area is less than or greater than the 
radius a of the junction between the regimes A and C. 
For 7) < a, equation (38) gives directly 


1 1 : 
FPR) +5 Ploy (0<r<»), | 


es Oy)? = | = p 1 log—— (SL Sa8)5 (47) 
5 on (@<7sR) 
For 7%) > a, the corresponding result is 
{b(@—P) (0<r<a), 
pip Wim 1 ple SE es eee (48) 
: pee sd Gee R) | 
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In each case the junction radius a must be-Such that equation (42) is satis- 
fied. The substitution of either form (47) or (48) for / into (42) gives a transcen- 
dental equation for a. The values of a/R were determined by graphical methods 
for the values 0, 0-2, 0-4, 0-6, 0-7,-0-8, 0-9 and 1:0 of 7/R. The results are 
shown graphically in Figure 3. The distances a and rg are equal when 7)/R is 
equal to 0-695 approximately. It can be seen from the figure that a > R/2 
and it follows from the discussion of inequality (41) that an associated collapse 
mode exists satisfying the condition (32) for a relative minimum volume 
design. 


0-75 
0-70 
a/R 
0-65 
0-60 
le} 0-2 0-4 0-6 0-8 me) 
r/R 
Figure 3 


Variation of a/R with 79/R for the solid built-in plate under circular loading. 


The minimum weight design for the case of uniform loading (7)/R = 1) is 
shown diagrammatically in Figure 2 (b) and may be contrasted with the sand- 
wich plate design [Figure 2 (a)] for the same conditions of loading and support. 

For a concentrated load P at the center, 7, = 0 and the thickness is given by 


: = log = (0 272%), 
— Fie = 49 
4 0 

IP y—a 

ee (a<r<R) 

UREA Y 


Here a has the value 0-74 R in order to satisfy equation (42). The volume of the 


plate is found to be : 
V = 0-407 (=) ~ R?. (50) 
0 


A solid plate of uniform thickness which is just at collapse under a central 


load P has the volume ae ; 
V, = 0-798 (=\" x R? (51) 
0 


for both the built-in and simply supported edge conditions. Comparison of 
(46) and (50) with (51) shows that the savings effected by designs (45) and (49) 
over the constant thickness design are 37°, and 49% respectively. 


ZAMP VIII/32 


498 Emin TuRAN ONAT, WALTER SCHUMANN and RicH. THORPE SHIELD ZAMP 


5. Built-in Sandwich Plate with Weight 


In [4], it was shown that, in principle, body forces may be taken into 
account in the optimum design of a structure. The design of a simply supported 
circular sandwich plate loaded by uniform pressure and its own weight was 
obtained as an illustration. Here as a further illustration we obtain the corres- 
ponding design for a built-in plate. 

The plate is assumed to be horizontal and loaded by a uniform pressure p 
on its upper surface. As in section 3, the core thickness H is constant and it 
can be supposed that the weight of the core is included in the pressure p. A 
minimum volume design is obtained if the plate is designed to collapse in a 
mode such that the rate of dissipation of energy due to plastic action less the 
rate at which the body force does work is constant in the face sheets (4). This 
condition can be written 


MnxzxtwNA 
a i a y w =const , (52) 


where y is the weight per unit volume of the material of the face sheets. 

We assume that regime A of Figure 1 applies for 0 <7 < a and regime C 
applies for a << 7 < R, where R is the radius of the plate. Condition (52) then 
becomes 


d2w 1 dw 
eta atk tee er ee | a 
53 
d? 
“iyi B72 w= o (Qa. 7 Se | 


where 6? = 2 y/o, H and « is a positive constant. Equations (53) are subject 
to the conditions w = dw/dr = 0 at the edge y = R of the plate and dw/dr is 
zero at the center of the plate. Also w and dw/dr are continuous at 7 = a. It is 
found that 


a [Jo(B 7) — Jo(B 2) sech B (R — a)] Os7 2 a5 
Wie | a Jo (8 a) [cosh B (R — 7) sech B (R — a) — sechB (R — a)] (54) 
(427k) 


and the radius a is given by 


tanh B (R — a) = ee (55) 


Equation (55) has a unique solution for a for given values of 6 and R, and the 


collapse mode (54) satisfies the appropriate inequalities on x and A in the inner 


and outer regions. 
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For equilibrium the moments satisfy the equation 


Q? Le aN 


a (y M) ip +rpt+2ayir=0. (56) 


In the region0 <r < a,M = N= Myandintheregiona <7 < R,M=-—M,, 
N = 0, where M, is given by (6). At the junction r = a, M is continuous and 
therefore zero and the shearing force Q is continuous. At the center Q is zero. 
It is found that 


ool i= B [cosh (ry—a)—1 (57) 
+ [4, + % tanh (R — a)] sinhp ( — a)} 
(a<r<R). 


As the parameter f increases, that is as the ratio of the weight of the 
material to its yield stress increases, the ratio a/R decreases, and the section 
y =a where the face sheets have zero thickness moves nearer to the center 
of the plate. 
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Zusammenfassung 


Diese Arbeit befasst sich mit dem Problem der Bemessung einer elastisch- 
plastischen Kreisplatte veranderlicher Dicke von minimalem Gewicht. Die Platte 
ist am Rand frei drehbar gelagert oder eingespannt und tragt eine rotations- 
symmetrische, sonst beliebig verteilte Last. Es wird das von DRUCKER und 
SHIELD entwickelte Verfahren zur Dimensionierung verwendet, welches ein abso- 
lutes Minimum fiir das Gewicht einer Sandwichplatte und ein relatives Minimum 
fiir dasjenige einer homogenen Platte liefert. 
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Engineering Analysis. A Survey of Numerical Procedures. Von 
STEPHEN H. CRANDALL (McGraw-Hill Book Company, New York, Toronto, 
London 1956). 417 S., 170 Fig.; €3.11s.6d. 

Die Aufgabe des vorliegenden Buches besteht darin, dem Ingenieur zu helfen, 
seine Probleme numerisch zu lésen. Die Fachrichtung spielt eine sekundare Rolle, 
da iiberall im Grunde genommen die gleichen mathematischen Betrachtungen 
anzustellen sind, wie Schaffung eines mathematischen Modells fiir eine physi- 
kalische Sachlage und die Reduktion des mathematischen Problems zu einem 
numerischen Verfahren. Vorwiegend wird der zweite Teil im Buch behandelt, 
doch gibt der Autor manche Beispiele fiir den ersten Schritt, in der Meinung, dass 
die Struktur der mathematischen Formulierung diskutiert werden muss zum 
Verstandnis der numerischen Behandlung. Drei Hauptprobleme sind es, die sich 
im allgemeinen dem Ingenieur darbieten: Gleichgewicht (eingeschwungener 
Zustand von Temperaturverteilungen, laminare Strémungen, elastische Span- 
nungszustande) ; Bestimmung von Eigenwerten (kritische Werte von bestimmten 
Parametern im Gleichgewichtszustand, kritische Frequenzen und Resonanzpro- 
bleme) und Ausbreitungsvorginge (Warme, Druckwellen in Fliissigkeiten und 
mechanische Spannungen und Verformungen in elastischen Systemen). Alle 
Probleme werden in Systemen mit beschranktem Freiheitsgrad (konzentrierte 
Parameter) sowie in kontinuierlichen behandelt. Auch nichtlineare Probleme, 
mit Ausnahme der Eigenwertbestimmung, kommen zur Sprache. Zahlreiche 
Ubungsaufgaben dienen der Festigung des Gebotenen. 

Dem rechnenden Ingenieur kann das Buch eine wertvolle Hilfe sein. 

H. Weber 


Wahrscheinlichkeitstheorie. Von H. RicHTER (Springer-Verlag, Berlin 1956 
[Grundlehren der Mathematischen Wissenschaften, Band 86]). 435 S., 14 Abb.; 
DM 66.—, Ganzleinen DM 69.60. 

Das hier vorliegende Lehrbuch ist das erste ausfiihrliche Werk in deutscher 
Sprache, das die neue Entwicklung der Wahrscheinlichkeitsrechnung beriick- 
sichtigt, und sein Erscheinen ist schon aus diesem Grunde sehr zu begriissen. 

Der Inhalt gliedert sich in die folgenden Abschnitte: Mathematische Grund- 
lagen, der Wahrscheinlichkeitsbegriff, die Elemente der Wahrscheinlichkeits- 
theorie, Elemente der Integrationstheorie, zufallige Grossen auf Wahrscheinlich- 
keitsfeldern, spezielle Wahrscheinlichkeitsverteilungen, die Konvergenz zufal- 
liger Grossen. Besonders ausfiihrlich werden die Grundlagen der Wahrscheinlich- 
keitstheorie behandelt, und zwar werden die Axiome ungefahr im Sinne von 
friiheren Ausfiihrungen des Verfassers (Math. Ann. 1953-1954) aufgestellt. Das 
vorgeschlagene Axiomensystem scheint besonders von Anwendungen in physi- 
kalischen Problemen beeinflusst zu sein. Die modernen Teilgebiete der Wahr- 
scheinlichkeitsrechnung, die stochastischen Prozesse, die Spieltheorie und die 
Theorie der statistischen Entscheidungsverfahren werden nicht behandelt. 

Vom Leser wird verlangt, dass er mit der mathematischen Denkweise gut 
vertraut ist, hingegen wird sehr wenig konkretes Wissen vorausgesetzt; insbeson- 
dere enthalt das Buch die notwendigen mengentheoretischen Voraussetzungen 


selbst. P. Schmid 
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Waffe und Wirkung bei der Fliegerabwehr. Von H. Branpti (Birk- 
hauser Verlag, Basel und Stuttgart 1956). 91 S., 24 Abb.; Fr. 26.—/ DM 26.-. 
Bei der vergleichenden Beurteilung verschiedener Flabgeschiitze ist es zufolge 
der grossen Zahl der auftretenden Parameter sehr schwierig, einen Gesamtiiber- 
blick zu erlangen. Der Verfasser des vorliegenden Buches versucht fiir den Fall 
des Folgeschiessens die theoretischen Grundlagen fiir einen solchen Vergleich 
herauszuarbeiten. In einem ersten Teil wird das sogenannte Prazisionsmass 
unter Beriicksichtigung der Geschossflugzeitstreuung untersucht und explizit 
dargestellt. Daraus kann dann der ‘Abschusskoeffizient, der bei Vernachlassigung 
wirtschaftlicher Uberlegungen taktisch allein massgebend ist, abgeleitet werden 
und zu einem Vergleich der Wirksamkeit verschiedener Waffen fiihren. Im 
zweiten Teil wird auf Grund gewisser Annahmen die Abhangigkeit der Treffwahr- 
scheinlichkeit von der Geschossflugzeit untersucht und ein gelegentlich verwen- 
deter Potenzansatz einer Kritik unterzogen. Schliesslich wird noch die beniitzte 
Zielfehlerdefinition gerechtfertigt. — Auch wenn jedenfalls noch verschiedene 
Hypothesen naher zu untersuchen bleiben, diirfte doch das Buch einen bedeut- 
samen Beitrag zur Klarung dieser ausserordentlich komplizierten Zusammen- 
hange liefern, Die iibersichtliche und klare Darstellung wird alle Fachleute inter- 

essieren, die mit schiesstechnischen Problemen zu tun haben. 
E. Roth-Desmeules 


Die Transformatoren. Von MILAN Vipmar, dritte, vollstandig neu bear- 
beitete Auflage (Birkhauser Verlag, Basel und Stuttgart 1956). 630 S., 321 Abb.; 
geb. Fr./DM 68.-, brosch. Fr./DM 64.-. 

Der Leser wird in acht Hauptabschnitten in lebendiger und klarer Weise 
griindlich mit den modernen Transformatoren bekannt gemacht. Es stellt ein in 
sich abgeschlossenes Werk dar, das sowohl Grundlagen und Eigenschaften wie 
Probleme beim Bau und Betrieb behandelt. Speziell ausfiihrlich werden die~ 
Hauptvertreter, die Drehstromtransformatoren, beschrieben, wahrend weniger 
ausgesprochene Transformatorenprobleme, wie Regulierapparate usw., nur er- 
wahnt bleiben. Hauptprobleme, die auch andere Elektromaschinenbauer inter- 
essieren, sind sehr breit behandelt, so die Warmeabfuhr aus dem arbeitenden 
Transformator, der Eisenkern, Leerlauf, Kurzschluss, Einschaltstromstoss und 
die Wicklung. Losgelést von allen Landes- und Firmennormen, wird vom physi- 
kalischen und wirtschaftlichen Standpunkt alles Vorhandene kritisch gegeneinan- 
der abgewogen und diskutiert sowie am Schlusse die Wege des Entwurfes anhand 
von ausgefiihrten Transformatoren illustriert. Das Buch wird sowohl dem Prak- 
tiker als dem Wissenschaftler und den Studenten gute Dienste leisten. 

G. Gubelmann und R. Weller 


Statistical Analysis of Stationary Time Series. Von U. GRENANDER und 
M. RosenBiattT (John Wiley and Sons, New York 1957). 300S., 8 Fig.; $11.-. 

Dieses bemerkenswerte Buch kann nur von solchen Lesern verstanden werden, 
welche gute Kenntnisse in der Wahrscheinlichkeitsrechnung, der mathematischen 
Statistik und der Masstheorie besitzen. Die Verfasser setzten sich das hohe Ziel, 
die statistische Analysis von Zeit-Reihen gemass modernen, exakten, mathema- 
tischen Methoden mit Hilfe der Masstheorie, Fourier-Analyse und Theorie der 
Matrizen darzustellen. Angesichts der grossen Bedeutung solcher Methoden in 
der modernen physikalischen und technischen Forschung sollten darin Tatige sich 
mit Hilfe dieses Buches iiber den betreffenden Problemkreis und seine mathe- 
matisch-statistische Behandlung orientieren. W. Sawer 
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Nonparametric Methods in Statistics. Von D. A. S. FRASER (John Wiley 
and Sons, New York 1957). 300-S., 17 Fig.; $8.50. 

In der modernen mathematischen Statistik werden neben klassischen Vertei- 
lungen, wie der Normalverteilung, definiert durch bestimmte Parameter, allge- 
meinere Verteilungen ohne Parameter betrachtet. Ebenso wurden in der stati- 
stischen Schatzungstheorie entsprechende Methoden geschaffen, die sich neben 
wahrscheinlichkeitstheoretischen Begriffen auf die Masstheorie stiitzen. Dieses 
Buch gibt einen guten Uberblick iiber den heutigen Stand dieser Theorie, zu 
welcher der Verfasser wichtige Beitrage beisteuerte. W. Saxer 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLiiccr. Bd. 20: Elektrische Leitungsphinomene II (Springer-Verlag, Berlin 
1957) 4Ote Sa 22ers. DMA TT ze 

Band 20 des Handbuchs behandelt elektrische Leitungsphanomene in Nicht- 
metallen und enthalt die vier Beitrage: O. MApELUNG, Halbleiter; A. B. Lip1arp, 
Ionic Conductivity; J. M. STEVELS, The Electrical Properties of Glass; und E. Dar- 
Mots, Electyochimie. 

Der Artikel von MADELUNG, der die Halfte des Bandes einnimmt, gehort zu 
den besten unter den neueren zusammenfassenden Arbeiten auf dem Gebiet der 
Halbleiterphysik. In einer bester Handbuchtradition entsprechenden Darstel- 
lung — es sei hier nur auf die Zusammenstellung der Koeffizienten der elektrischen, 
magnetischen und thermischen Effekte in Abschnitt C hingewiesen — wird ein 
Uberblick iiber das weitschichtige Gebiet gegeben. Dass dabei die Diskussion 
spezieller Halbleiter etwas zu kurz kommt, diirfte zum Teil auf Platzmangel, zum 
Teil aber auch auf die gegenwartig noch sehr liickenhaften experimentellen Ergeb- 
nisse zurtickzufiihren sein. 

Auch die zweite, von Lipiarp stammende Arbeit erfiillt in vollem Umfang die 
Anspriiche, welche an einen Handbuchartikel gestellt werden. In kurzer und 
tibersichtlicher Form werden auf 100 Seiten alle wesentlichen Gebiete der Ionen- 
leitung behandelt. Besondere Erwahnung verdient eine Zusammenstellung neuer 
und neuester Literatur, die, nach Substanzen getrennt, vom Fachmann sehr 
begriisst werden diirfte. 

Entsprechend dem heutigen Stand der Forschung enthalt der Beitrag von 
STEVELS iiber die elektrischen Eigenschaften von Glas vorwiegend experimentelle 
Ergebnisse. 

Man fragt sich, ob die Arbeit von Darmots iiber Elektrochemie im 20. Band 
des Handbuches an der richtigen Stelle steht. (Der Artikel von GarRLicxk, Bd. 19, 
iiber Photoleitung hatte unseres Erachtens besser in diesen Zusammenhang 
gepasst.) Im tbrigen lasst diese Arbeit in bezug auf Ubersichtlichkeit sowohl als 
auch in der Berticksichtigung neuester Ergebnisse zu wiinschen iibrig. E. Mooser 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLucce. Band 32: Strukturforschung (Springer-Verlag, Berlin 1957). 663 S., 
373 Fig.; DM 144.-. 

Es ist erstmalig, dass alle Arten von Strukturuntersuchungen mit Roéntgen-, 
Elektronen- und Neutronenstrahlen in einem einzigen Bande dargestellt werden. 
Nicht nur das Bestimmen von Kristallstrukturen, sondern auch Untersuchungen 
von Fliissigkeiten, amorphen und makromolekularen Stoffen werden darin be- 
handelt. Der schén ausgestattete und reich illustrierte Band enthalt wertvolle 


Literaturhinweise auf neueste Arbeiten und ist als Hand- und Lehrbuch deshalb 
sehr willkommen. 
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Der einfiihrende Artikel von A. GuInigrr“und G. v. ELLER, Les méthodes 
expérimentales des déterminations de structures cristallines pay vayons X, beschreibt 
ausser den experimentellen Methoden die einzelnen Arbeitsschritte einer voll- 
standigen Kristallanalyse besonders klar. Die theoretische Begriindung dazu, 
insbesondere die Theorie der Réntgénstreuung und -reflexion, wird im Artikel 
von J. Bouman, Theoretical Principles of Structural Research by X-Rays, gegeben. 
G. FouRNET zeigt in Etude de la structure des fluides et des substances amorphes au 
moyen de la diffraction des rayons X durch genaue theoretische Untersuchung, 
dass die alteren Deutungen von Fliissigkeitsinterferenzen nur teilweise richtig 
sind. Die Kleinwinkel-Réntgenstreuung findet interessante Anwendungen in der 
Virus- und Eiweissforschung, wie aus der Arbeit von W. W. Beeman, P. KagEs- 
BERG, J. W. ANDEREGG und M. B. Wess, Size of Particles and Lattice Defects, 
hervorgeht. Im zweiten Teil wird gezeigt, welche Aussagen iiber Gitterfehler 
durch réntgenographische Untersuchungen gewonnen werden kénnen.-In den 
Beitragen von H. RAETHER iiber Elektroneninterferenzen und von G. R. RInGo 
liber Neutron Diffraction and Interference werden besonders die Unterschiede 
und die Vorziige gegeniiber den R6ntgenverfahren fiir bestimmte Strukturpro- 
bleme hervorgehoben. H.. Grdnicher 


Physikalisch-statistische Regeln als Grundlagen fiir Wetter- und 
Witterungsvorhersagen. Von FRANz Baur. 1. Band (Akademische Verlags- 
gesellschaft mbH., Frankfurt am Main 1956). 139 S., 30 Abb., 2 farbige Tafeln; 
geb. DM 46.-. 

In diesem Werk, dessen erster Band jetzt vorlegt, unternimmt es der Verfas- 
ser, auf Grund jahrzehntelanger Forschungen und anhand eines umfangreichen 
Materials nachzuweisen, «dass die physikalisch durchdachte, wahrscheinlichkeits- 
theoretisch kontrollierte Statistik ganz allgemein die den Wettervorgangen am 
besten angepasste Form der Erfahrungsforschung ist und dass mit ihr auch fiir 
die kurz- und mittelfristige Wettervoraussage bessere und tiberzeugendere Er- 
fahrungsgrundlagen gewonnen werden konnen als mit den bisher fast ausschliess- 
lich gebrauchten Beispielen vereinzelter Wetterlagen». Dieser Nachweis ist dem 
Verfasser bis zu einem hohen Grade tatsachlich auch gelungen. Es konnte eine 
ganze Reihe zutreffender Regeln iiber Zusammenhange zeitlich sich folgender 
Witterungserscheinungen abgeleitet werden, die auch fiir weitere theoretische 
und empirische Untersuchungen auf dem Gebiet der Witterungsvorhersage als 
Ausgangspunkt dienen k6énnen. 

Im ersten Abschnitt werden die statistischen Grundlagen fiir die Vovaussagen 
im tdglichen Wetterdienst auf das sogenannte Azorenproblem angewandt, das 
heisst auf die Beziehung zwischen Luftdruckanderungen auf den Azoren und 
nachfolgenden Druckanderungen in Mitteleuropa mit den sie begleitenden Wet- 
tererscheinungen. Vorhersagen, die nicht nur fiir lange Zeitabschnitie, sondern 
auch fiir Grosswetterlagen gelten sollen, kénnen nur bei sehr genauer Kenntnis der 
Verhiltnisse der allgemeinen atmospharischen Zirkulation gemacht werden. Der 
zweite Abschnitt befasst sich daher mit theoretischen Betrachtungen iiber die 
Schwankungen der allgemeinen Zirkulation als Richtschnur fir die statistische 
Grosswetterforschung, wahrend im dritten Abschnitt die Beobachtungstatsachen 
der taglichen atmospharischen Zirkulation auf der Nordhalbkugel in den Jahren 
1949-1951 dargestellt und gedeutet werden. Schliesslich wird im vierten Ab- 
schnitt das ebenso weitschichtige wie schwierige Problem der Zusammenhange 
des Grosswetters mit solaren Vorgangen behandelt, indem Regeln tiber die Zu- 
ordnung des Grosswetters zum Sonnenfleckenzyklus gegeben werden. 
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Der Synoptiker, dessen Voraussagen auf sorgfaltigen Untersuchungen ein- 
zelney Wetterlagen, auf einem guten Gedachtnis und nicht zuletzt auf dem Finger- 
spitzengefiihl aufgebaut sind, kann von den statistischen Methoden der Vorher- 
sage nur profitieren. Das aus reicher Erfahrung entstandene Werk von F. Baur 
gehort ohne Zweifel zu den bemerkenswerten Erscheinungen der meteorologischen 
Literatur. J.C. Thams 


Fortschritte in der meteorologischen Forschung seit 1900. Von BERn- 
HARD Nets (Akademische Verlagsgesellschaft mbH, Frankfurt am Main 1956). 
238 S., 21 Abb., 1 Titelbild und 3 Portrattafeln; geb. DM 28.-. 

Dieses Werk ist aus Vorlesungen an der Freien Universitat Berlin hervor- 
gegangen. Es ist weder ein Lehrbuch noch ein Geschichtswerk. Aus dem vor- 
wiegend deutschsprachigen Schrifttum wahlt es die Arbeiten aus, die wesentlich 
dazu beigetragen haben, die Meteorologie zu einer exakten Naturwissenschaft 
auszugestalten. 

Das Werk gliedert sich in vier Abschnitte. Im ersten wird auf die Leitideen der 
meteorologischen Forschung eingegangen; im zweiten wird das aerologische 
Beobachtungsgut gedeutet. Der dritte Abschnitt behandelt das Wetter als eine 
Folge der Energieumsatze in der Atmosphare, wahrend der vierte Abschnitt der 
praktischen Meteorologie gewidmet ist. Als Ausgangspunkt wird der Anfang des 
20. Jahrhunderts genommen, den der Verfasser als einen Wendepunkt in physi- 
kalischer und meteorologischer Hinsicht betrachtet, da erst von 1900 an die 
Arbeiten geleistet worden seien, welche die Aufstellung einer Theorie des Wetters 
ermoglichten. 

Ohne Zweifel entspricht eine solche kritisch-genetische Schau einem Bediirf- 
nis, und es ist gewiss verdienstvoll, wenn der Verfasser versucht, die treibenden ~ 
Motive in der Entwicklung der wissenschaftlichen Meteorologie aufzuzeigen. Inso- 
fern kann das Werk wirklich ein Helfer beim Meteorologiestudium sein. Es hatte 
jedoch erheblich gewonnen durch Beriicksichtigung der gesamten einschlagigen 
Literatur. J.C. Thams 


Acht- und neunstellige Tabellen zu den elliptischen Funktionen. Von 
M. ScHULER und H. GEBELEIN (Springer-Verlag, Berlin 1955). 296 S., 11 Abb.; 
DM 58.—. 

Funfstellige Tabellen zu den elliptischen Funktionen. Von M. ScHULER 
und H. GEBELEIN (Springer-Verlag, Berlin 1955). 114 S., 11 Abb.; DM 29.60. 

Beide Werke sind mit vollstandigem deutschem und englischem Text verse- 
hen und vor allem auf die Méglichkeit einer bequemen Interpolation ausgerichtet. 
Inhalt: 

Jacobische elliptische Funktionen, laufend nach z und q. 

Hilfsfunktionen fiir eine einfache Berechnung der #-Funktionen. 

Tafeln fiir die Umrechnung zwischen dem Legendreschen Modul @ und dem ~* 

Jacobischen Parameter gq. 
Im fiinfstelligen Werk ferner: Koeffizienten fiir die Everett-Interpolation. 


Tables of Weber Parabolic Cylinder Functions. Von J.C. P. MILLER 
(Her Majesty’s Stationery Office, London 1955) 0233) S3n03Ks 

Tabelliert sind die Lésungen W(a, x) und W(a, —x) der Differentialgleichung 

y" + (4#*/4 — a) y= 0 mit reduzierten Ableitungen, fiir a= —10(1)10 und 

« = 0 (0,1) 10, sowie verschiedene Hilfsfunktionen. Ausfiihrliche mathematische 

Einleitung, welche unter anderem die Querverbindung zu anderen Funktionen 

herstellt. P. Laéuchli 
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